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^ i ABSTRACT 

^ ; The AdS/CFT correspondence allows us to map a dynamical cosmology to a dual quantum 
field theory living on the boundary of spacetime. Specifically, we study a five-dimensional 
model cosmology in type IIB supergravity, where the dual theory is an unstable deformation 
of Af = 4 supersymmetric SU(N) gauge theory on H x S 3 . A one-loop computation shows 
that the coupling governing the instability is asymptotically free, so quantum corrections 
cannot turn the potential around. The big crunch singularity in the bulk occurs when a 
boundary scalar field runs to infinity, in finite time. Consistent quantum evolution requires 
that we impose boundary conditions at infinite scalar field, i.e. a self-adjoint extension of 
the system. We find that quantum spreading of the homogeneous mode of the boundary 
scalar leads to a natural UV cutoff in particle production as the wavefunction for the ho- 
mogeneous mode bounces back from infinity. Translating back to the bulk theory, we find 
that a quantum transition from a big crunch to a big bang is the most probable outcome 
of cosmological evolution, for a specific parameter range. Intriguingly, the instability and 
approximate scale-invariance of the boundary theory lead to the generation of an approx- 
imately scale-invariant spectrum of stress-energy perturbations on the boundary, whose 
amplitude is naturally small. We comment on qualitative differences with holographic de- 
scriptions of large black holes, on four- dimensional generalizations and on implications for 
cosmological perturbations. 
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1 Introduction and Summary 

No theoretical framework for cosmology can claim to be complete until it resolves the 
cosmological singularity, most likely through quantum gravitational effects. We do not 
know, for example, whether the singularity represents the beginning of the universe and, if 
not, what happens to the thermodynamic arrow of time there. The no-boundary proposal 
[1], for example, predicts the arrow of time reverses. Another, perhaps simpler possibility is 
that the physics responsible for resolving the singularity also describes dynamical evolution 
across it [2-4]. 

Whichever of these options is correct has profound implications for cosmology. If the 
singularity was the beginning, the horizon and flatness problems, and the problem of the 
origin of the observed density variations, seem to demand an early epoch of inflation [5]. 
But why the universe emerged dominated by inflationary energy, yet with a tiny, fine- 
tuned cosmological constant is hard to explain. Whereas if the singularity was not the 
beginning and there was a preceding epoch of slow cosmological contraction (in Einstein 
frame) before it, non-inflationary solutions of the classic cosmological problems are possible 
[3,6,7]. Furthermore, in a universe where repeated cycles of evolution occur, technically 
natural mechanisms whereby the cosmological constant can relax to tiny values become 
viable [8]. 

The AdS/CFT correspondence [9] has emerged as an extremely powerful tool for un- 
derstanding quantum gravity, providing a non-perturbative definition of string theory in 
asymptotically anti-de Sitter (AdS) spacetimes, in terms of conformal field theories (CFTs) 
on their conformal boundaries. In this paper, we shall use the AdS / CFT correspondence to 
provide a quantum resolution of cosmological singularities. We study Maldacena's original 
AdS 5 x S 5 example in type IIB supergravity, but with generalized boundary conditions on 
some of the negative mass squared scalars (saturating the Breitenlohner-Freedman bound). 
These generalized boundary conditions allow smooth, asymptotically AdS initial data to 
evolve into a big crunch singularity, namely a spacelike singularity that reaches the space- 
time boundary in finite time [10, 11]. 

The dual description of these 'AdS cosmologies" involves field theories with scalar po- 
tentials which are unbounded below and which drive certain boundary scalars to infinity in 
finite timely The AdS/CFT duality, therefore, relates the problem of resolving cosmological 
singularities to that of understanding the dynamics of quantum field theories of this type. 

In our example, the dual field theory is a deformation of M = 4 Super- Yang-Mills (SYM) 
theory on Kx S 3 by an unbounded double trace potential —f0 2 /2, with O a trace operator 
quadratic in the adjoint Higgs scalars [12] . The dual field theory has several key properties 
that allow us to analyze its dynamics quantitatively: 

1. The principal advantage of this five-dimensional setup, compared with similar four- 
dimensional cosmologies studied previously, is that the undeformed dual boundary 
theory is well-understood. Furthermore, the deformation is renormalizable and a 

1 There is strong evidence that supergravity solutions with big crunch singularities are always related to 
dual field theories of this form. This is related to cosmic censorship in the bulk [10]. 
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one-loop computation shows that, because of the "wrong" sign of the deformation, 
the coupling / that governs the instability is asymptotically free [13], allowing per- 
turbative field theory computations in the regime near the singularity, i.e., for large 
O, at least at small 't Hooft coupling. In particular, we can show that quantum 
corrections do not turn the potential around, so that it is really unbounded below. 

2. As the dual theory evolves towards the singularity, the field evolution becomes ul- 
tralocal on any fixed length scale, meaning that spatial gradients become dynamically 
unimportant. This means that one effectively has an infinite set of decoupled quantum 
mechanical systems, one at each spatial point, in the regime near the singularity. 

3. Since the conformal boundary has finite spatial volume, the unstable, homogeneous 
background mode evolves quantum mechanically. The quantum mechanical spread of 
its wave function will turn out to be crucial for the suppression of particle creation as 
the wavefunction for the homogeneous mode rolls down the potential, bounces back 
from infinity, and rolls back up towards its starting value. 

4. As the singularity is approached, the semiclassical approximation becomes increas- 
ingly accurate, both for the homogeneous background and for the fluctuations. This 
allows us to study the full quantum dynamics with some analytical precision. 

In the unstable dual theory, even if we start with the homogeneous background mode 
described by a localized wave packet, the wavefunction spreads to infinite scalar field in an 
arbitrarily short time. Unless suitable boundary conditions are imposed, probability will be 
lost at infinity. Therefore, in order to study quantum evolution in this theory at all, one has 
no choice but to impose unitary boundary conditions at large field values, i.e., boundary 
conditions which restrict the Hilbert space to a subspace on which the Hamiltonian is self- 
adjoint. In quantum mechanics, such a restriction is known as a "self-adjoint extension" 
of the original ill-defined theory [14,15]. Because the field theory evolution is ultralocal 
near the singularity, we can attempt to extend this construction to the full boundary field 
theory by imposing the same self-adjoint extension at each spatial point. 

The fact that the semiclassical approximation becomes more and more accurate near 
the singularity allows us to implement the self-adjoint extension using classical solutions 
and the method of images. The relevant solutions turn out to be generically complex near 
the singularity as a consequence of the quantum spread in the homogeneous mode. To lead- 
ing order, therefore, the inhomogeneous modes evolve in a complex classical homogeneous 
background. The complexity of the background turns out to be essential to the resolution of 
the singularity, providing an ultraviolet cutoff on quantum particle creation when the wave 
packet rolls down the potential and bounces back. We shall find that, for a certain range 
of parameters, backreaction of created particles on the homogeneous mode is negligible for 
values of the final Schrodinger wavef unction's argument carrying significant probability. 

Our procedure specifies a consistent rule for unitary quantum evolution in the boundary 
theory, in the presence of bulk cosmological singularities. Evolving the wave packet forward 
in time, the most probable outcome is that the homogeneous background returns close to 
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its original starting point, with a finite density of particles created at the bounce. Therefore 
we predict a transition from a big crunch to a big bang in the bulk as the most probable 
outcome of cosmological evolution, for a certain range of parameters. 

If the 't Hooft coupling in the boundary super- Yang-Mills theory is small, we have good 
control over the field theory. However, for small 't Hooft coupling the bulk is in a stringy 
regime. In fact, we shall see that the form of the quantum effective potential remains valid 
at large 't Hooft coupling, so at least the key dynamical feature of an unbounded negative 
potential is shared by both regimes. Further analysis is required, though, to fully treat the 
boundary theory in the regime where the bulk is well- described by supergravity. 

Let us now briefly review the status of other approaches to cosmological singularity 
resolution. The simplest string theory models with cosmological singularities are time- 
dependent orbifolds. Here, however, perturbative string theory tends to break down, with 
large gravitational backreaction leading to divergences in perturbative string scattering 
amplitudes [16]. It is natural to hope that a condensation of winding modes can improve 
the situation, but in general this is still unclear [17]. However, in a specific model it has 
been argued that the singularity is replaced by a winding tachyon condensate phase, and 
that the system can be analyzed within perturbative string theory [18]. 

The ekpyrotic model [3] involves similar time-dependent orbifold solutions to M-theory. 
Near the singularity, the only light modes in the theory are winding M2-branes [19,20], and 
the effective string coupling (in IIA or heterotic frame) vanishes. A study of the classical 
M2-brane dynamics suggests that the a' expansion must be replaced by an expansion in 
1/a' [19], but fully quantum mechanical calculations have not yet been possible. 

Other attempts to describe spacelike singularities using a non-perturbative dual de- 
scription have also been made. Models with two spacetime dimensions have been described 
using c = 1 matrix models [21]. Higher-dimensional models with light-like singularities 
have been studied in the framework of matrix theory [22], which has led to the suggestion 
that spacetime may be replaced by non-commuting matrices near a singularity [23]. The 
quantitative study of the dynamics of this regime is still work in progress, however. 

In earlier work, the AdS/CFT correspondence has also been used to study the singularity 
inside black holes, which is closely analogous to a cosmological singularity [24]. Although 
some progress in this direction has been made, the fact that the singularity is hidden 
behind an event horizon clearly complicates the problem. The CFT evolution is dual to 
bulk evolution in Schwarschild time, so the CFT never directly "sees" the singularity. This 
should be contrasted with the model discussed in the present paper, in which the bulk 
singularity reaches the boundary in finite time and is thus directly visible in the boundary 
theory. Other AdS / CFT models of cosmological singularities include [25] . 

At first sight our findings seem to indicate that the quantum dynamics of cosmological 
singularities is rather different from the thermalization process that describes the formation 
of large AdS black holes in the dual theory. To compare the two situations, we show the 
Penrose diagram of gravitational collapse to a large black hole in anti-de Sitter space in 
Figure [TJ This is identical to the Penrose diagram of an anti-de Sitter cosmology. There- 
fore from this point of view, there appears to be little distinction in the classical bulk 
theory between singularities inside black holes and cosmological singularities extending to 
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big bang 



Figure 1: The Penrose diagram of a large black hole in anti-de Sitter space is identical to 
that of an anti-de Sitter cosmology. In a black hole spacetime, however, time at infinity 
continues forever whereas in an AdS cosmology the singularity hits the boundary in finite 
time. 

the boundary in finite time. However the dual (quantum) description of both types of 
singularities appears to be qualitatively different: as we discussed above, AdS cosmologies 
are described in terms of unstable conformal field theories with steep potentials V that are 
unbounded below (Figure [21 left), whereas black holes are interpreted as thermal states 
requiring a vacuum state in the dual theory (Figure [21 right). 

If one "regularizes" the unbounded potential in Figure [2] (left), for example by adding 
higher order terms to obtain a potential like the one shown in Figure [2] (right), one finds 
this changes the evolution in the bulk near the big crunch, in the upper corners of the 
Penrose diagram. This is because the regularization affects the bulk boundary conditions. 
In particular one finds this causes the cosmological singularity to turn into a large (stable) 
black hole with scalar hair [11], where the bulk scalar field turned on is dual to the operator 
O in the boundary theory. This is a new type of black hole which does not exist (in a 
stable form) for the original bulk boundary conditions. It has a natural interpretation in 
the dual theory as an oscillatory excitation about the global negative minimum of V [11] 
(whereas the usual Schwarschild-AdS black holes correspond to thermal states around the 
standard vacuum at (O) = 0). This new black hole with scalar hair is the natural end 
state of evolution in the bulk corresponding to a wave packet rolling down a regularized 
potential. 

Considering a series of dual theories where the global minimum is taken to be more and 
more negative, one finds that the horizon size of the black holes with scalar hair - keeping the 
mass constant - increases. In the limit where the global minimum goes to minus infinity the 
hairy black holes become infinitely large and we recover the original cosmological solutions. 
While for the black hole case (with bounded potentials) it is plausible that the system will 
eventually thermalize, it is unclear what the very late time behaviour is in the cosmological 
case (with unbounded potentials). But the two cases may be similar for intermediate 
times. Indeed one expects that even for boundary conditions corresponding to steep but 
bounded potentials a wave packet rolling down the potential will bounce back a number of 
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Figure 2: The dual description of AdS cosmologies involves a wave packet rolling down 
an unstable direction of the field theory potential (left). At first sight this is qualitatively 
different from the formation of large black holes in AdS, which is described as a thermal- 
ization process in a dual theory that has a ground state (right). We argue however that for 
cosmology the dynamics at intermediate times is relevant, and that this may well be qual- 
itatively similar in both theories, at least for potentials that are sufficiently steep around 
the global minimum. 



times before the system thermalizes and settles down. These bounces could presumably be 
interpreted as "regularized" versions of the big crunch/big bang transitions we will study. 

We emphasize, however, that although the dual description in terms of a bounded V 
describing the formation of hairy black holes may share some qualitative features with the 
cosmologies at intermediate times, it does not provide us with a precise and consistent 
model of singularities. This is because the higher order terms needed to regularize the 
potential generally lead to a non-renormalizable boundary theory. In contrast, in this 
paper we focus on a class of five-dimensional cosmologies for which the dual field theory 
is a renormalizable deformation of Af = 4 SYM on K x S 3 by an unbounded double trace 
potential —f0 2 /2. As we have explained, in this model we are able to develop a consistent 
and precise description of at least some quantities of interest, described by unitary quantum 
evolution in the boundary theory in the presence of bulk cosmological singularities. 

Although the prime focus of this work is on understanding whether passage across a 
cosmological singularity is possible, there is another, perhaps even more remarkable con- 
sequence of the AdS/CFT dual description of cosmology. Since the boundary theory is 
unstable, and undergoes a bounce, a spectrum of fluctuations is generated by the ekpy- 
rotic mechanism [26]. As the unstable scalar rolls down the potential, the inhomogeneous 
modes acquire an increasingly negative mass squared, causing modes of successively higher 
momentum to enter a "growing mode" solution. By this process, a spectrum of spatially 
inhomogeneous excitations is generated across the bounce. The boundary theory is classi- 
cally scale-invariant. Quantum effects break this symmetry, but only slightly because the 
coupling / is asymptotically free. As a consequence of the approximate scale-invariance of 
the dual theory, the stress-energy fluctuations turn out to be approximately scale-invariant. 
Furthermore, their amplitude is suppressed by positive powers of the small parameter /. 
While we still have to translate these boundary fluctuations into the bulk, it is intriguing 
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that the AdS/CFT duality seems to automatically yield a physical mechanism for producing 
scale-invariant, small-amplitude cosmo logical perturbations. 

The outline of this paper is as follows. In Section 2, we introduce the bulk cosmologies 
of interest. We review how modifications of the standard AdS boundary conditions allow 
smooth initial data to evolve into a big crunch singularity, and focus on a specific example 
for which the dual field theory analysis will be tractable. We also show that the bulk 
cosmological solution to IIB supergravity is, after a duality taking us to type IIA frame, 
qualitatively similar to that describing compactified Milne (or colliding orbifold planes) 
in M-theory. In Section 3, we discuss the dual field theory, an unstable double trace 
deformation of M = 4 SYM theory, deriving the effective potential both at weak and 
strong 't Hooft coupling. Section 4 discusses self-adjoint extensions in quantum mechanics, 
as well as the ultralocality that will allow us to extend this idea to quantum field theory. 
In Section 5, we discuss the quantum evolution of the homogeneous background, which 
exhibits a quantum spread because the field theory lives on a finite volume space. We 
develop a method for implementing self-adjoint extensions in the semiclassical expansion, 
using complex classical solutions and the method of images. We also explain how this 
method extends to potentials with branch points, such as the quantum effective potential 
of interest in this paper. Section 6 focuses on the inhomogeneous modes, in particular on the 
question whether abundant particle creation may prevent the scalar field from running back 
up the potential after the bounce. We find that for most of the range of the Schrodinger 
wavefunction's argument (pj carrying significant probability, the quantum spread of the 
homogeneous background provides an ultraviolet cutoff on the wavelength of produced 
particles; provided the spread is sufficiently large, backreaction on the homogeneous mode is 
suppressed. In Section 7, we calculate the two-point correlator of stress-energy fluctuations 
in the boundary theory, which in principle determines the spectrum of classical metric 
perturbations in the bulk. Appendix A contains more details on the bulk theory. In 
Appendices B and C, we discuss a number of technical details related to the field theory 
effective potential. Appendix D contains details on complex classical solutions with and 
without a cutoff in the field theory. 

A brief and less technical overview of this work can be found in [27]; a preliminary 
report was presented in [28]. 

2 Anti-de Sitter Cosmology 
2.1 Setup 

Our starting point is Af = 8 gauged supergravity in five dimensions [29-31], which is 
thought to be a consistent truncation of ten-dimensional type IIB supergravity on S* 5 . The 
spectrum of this compactification involves 42 scalars parameterizing the coset E 6 /^/USp(8). 
We concentrate on the subset of scalars that parameterizes the coset SX(6, R)/SO(6). From 
the higher-dimensional viewpoint, these arise from different quadrupole distortions of S 5 . 
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The relevant part of the action involves five scalars q.{ and takes the form [32] 



S 



(2.1) 



where we have chosen units in which the coefficient of the Ricci scalar is |, i.e., the 5d 
Planck mass is unity. The potential for the scalars «j is given in terms of a superpotential 
W(«j) via 
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The potential reaches a negative local maximum when all the scalar fields ccj vanish. 
This is the maximally supersymmetric AdS state, corresponding to the unperturbed S 5 in 
the type IIB theory. At linear order around the AdS solution, the five scalars each obey 
the free wave equation with a mass that saturates the Breitenlohner-Freedman (BF) bound 
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Nonperturbatively, the fields couple to each other and it is generally not consistent to set 
only some of them to zero. However it is possible to truncate this theory further [30] to 
gravity coupled to a single SO (5)-invariant scalar ip by setting^ = (p/\/30 for i = 1, .., 5 
and P G = —5p/\/30. The action (12.11) then reduces to 



S 



-R 



:(V*>) S 



4R 2 Ads 



(I5e 2w + We 



3 -K>7</> 



(2.6) 



with 7 = y/2/15. 



2 There are several inequivalent ways in which this theory can be further truncated to a single scalar as 
only matter field [32]. An alternative option that was studied in [10,11] is to take (3i = tp/V3 for i = 1, ..,4 
and [3$ = /?6 = — 2</j/v3, which corresponds to taking only 015 7^ 0. This choice preserves an 5*0(4) x 5*0(2) 
symmetry. 
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2.2 Boundary Conditions 



We will work mainly in global coordinates in which the AdS§ metric takes the form 

dsl = R\ dS (-(1 + r 2 )dt 2 + ^ + r 2 ^ 3 ) • (2.7) 
In all asymptotically AdS solutions, the scalar ip decays at large radius as 

<p{r) = — + ^, (2.8) 

where a and 8 generally depend on the other coordinates. 

For the dynamics of the theory to be well-defined it is necessary to specify boundary 
conditions at r = oo on the fields. This amounts to specifying a relation between a and 8 
in (12.81) . For example, one can take a = 0, leaving 8 totally unspecified. This is the usual 
boundary condition, which preserves the full AdS symmetry group and which has empty 
AdS as its stable ground state [34,35]. Alternatively, one can adopt boundary conditions 
of the form 

— W 

where W(8) is an essentially arbitrary real smooth functionjf) Theories of this type have 
been called designer gravity theories [37], since their dynamical properties depend signifi- 
cantly on the choice of W. We will see that under the AdS/CFT duality, this function W 
appears as an additional potential term in the action of the dual field theory. 

Boundary conditions of the form (12. 9p generically break some of the asymptotic AdS 
symmetries, but they are invariant under global time translations. The conserved energy 
associated with this is well-defined and finite [38,39], but its expression depends on the 
function W . 

This can be seen as follows. When W is nonzero, the scalar field falls off more slowly 
than usual. This backreacts on the asymptotic behavior of the g rfl metric components, 
which causes the usual gravitational surface term of the Hamiltonian to diverge. This 
divergence is exactly canceled, however, by an additional scalar contribution to the surface 
terms. The total charge can therefore be integrated (provided one has specified a functional 
relation between a and 8). Hence one arrives at a finite expression for the conserved mass, 
which generally contains an explicit finite contribution from the scalar field which depends 
on W. Whether or not the energy admits a positive mass theorem^] however, depends on 
the choice of the function W. 

Below we will be interested in solutions with the following scalar field boundary condi- 
tions 

<*f = fP, (2.10) 
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In general, boundary conditions specified by an arbitrary function W can be imposed in anti-de Sitter 
gravity coupled to a tachyonic scalar with mass m 2 in the range m 2 BF < m 2 < m 2 BF + R^ds (see e.g. [36]). 
4 See [40] for recent work on the stability of theories of this type. 
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Figure 3: Regular initial data (p(r) that evolve to a big crunch singularity for boundary 
conditions with / = 0.1. 



where / is an arbitrary constant. The corresponding asymptotic form of the g rr component 
of the metric is given by 
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The conserved mass of spherically symmetric configurations with these boundary conditions 
reads 

"3 „o A 1 



M = 2n 2 R\ dS 



-M + f3 2 



(2.12) 



2 \ 2 

where Mq is the coefficient of the 0(l/r 6 ) correction to the g rr component of the metric. 



2.3 AdS Cosmologies 

We now construct a class of asymptotically AdS big bang/big crunch cosmologies that are 
solutions of (12. 6p with boundary conditions (12 . 10[) on the scalar field, with / > 0. This is 
a straightforward generalization of the four- dimensional cosmologies discussed in [10,11]. 

A particularly simple example of an open FLRW cosmology can be found from the 
evolution of an initial scalar field profile ip(r) obtained from an 0(5)-invariant Euclidean 
instantorH- Indeed, in Appendix A we show that all boundary conditions (I2.10p . for / > 0, 
admit precisely one such instanton solutionis The slice through the instanton obtained 
by restricting to the equator of the four sphere defines time symmetric initial data for a 
Lorentzian solution with mass M = —ir 2 R 2 idS f 2 f3 2 /4. The instanton, therefore, specifies 
negative mass initial data in this theoryQ In Figure [3] we show the initial scalar field profile 
obtained in this way for / = 0.1. 



5 See e.g. [10,11,41] for a discussion of similar four-dimensional cosmologies. 

6 Instantons also exist for negative / provided / < —0(1). However, as we explain below, we do not 
expect initial data obtained from slicing these instantons across the four sphere to evolve to a big crunch. 

7 As mentioned earlier, the mass of initial data obtained from instantons depends on the asymptotic 
behavior of the fields. For the AdS-invariant boundary conditions discussed in Appendix A (see eq. IA.5I) 
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Figure 4: Anti-de Sitter cosmology. To predict what happens at the singularities one must 
turn to the dual field theory description. 

Analytic continuation of the Euclidean geometry yields a Lorentzian solution that de- 
scribes the evolution of these initial data under AdS-invariant boundary conditions [42]. 
The origin of the Euclidean instanton then becomes the lightcone emanating from the ori- 
gin of the Lorentzian solution. Outside the lightcone, the scalar field is constant along 
four-dimensional de Sitter slices of AdS and the scalar remains bounded in this region. On 
the light cone we have tp = ip{r = 0) and dt<p = (since ip ir = at the origin in the 
instanton). Inside the lightcone, the SO(4, 1) symmetry ensures that the solution evolves 
like an open FLRW universe, 

ds 2 = -dt 2 + a 2 {t)da A , (2.13) 

where da4 is the metric on the four- dimensional unit hyperboloid. Under time evolution, 
ip rolls down the negative potential. This causes the scale factor a(t) to vanish in finite 
time, producing a singularity that extends to the boundary of AdS in finite global time. A 
coordinate transformation in the asymptotic region outside the light cone between the usual 
static coordinates (12. Tp for AdS^ and the 5*0(4, 1) invariant coordinates (see Appendix A) 
shows that (3 ~ (3{t = 0)/(cost) 2 . Hence f3 is now time dependent and blows up as t — > 7r/2, 
when the singularity hits the boundary. 

The boundary conditions (12.101) of interest here break conformal invariance weakly when 
the logarithm is large, and one cannot obtain the Lorentzian solution by analytic continua- 
tion from an 0(5)-invariant instanton. Instead one must evolve the initial data numerically. 
However, when the logarithm is large, our boundary conditions are nearly AdS-invariant 
except when a becomes large, and causality then restricts its effect. In particular, since 

one finds the instanton initial data have exactly zero mass, in line with their interpretation as the solution 
AdS 5 decays into [10]. 
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the evolution of the instanton data for AdS-invariant boundary conditions has trapped sur- 
faces, a singularity will still form in the central region and the effect of the modification 
of the boundary conditions on the evolution will only be appreciable in the corners of the 
conformal diagram where the singularity hits the boundary at infinity. Hence it is reason- 
able to expect that a singularity must form under evolution with aj boundary conditions 
when / is small. This restricts us to positive / since for negative / we find instantons only 
for / < —O(l). It would now appear possible, however, for the singularity to be enclosed 
inside a large black hole instead of extending out to infinity. Since the instanton initial 
data have slightly negative mass, MRj^ dS ~ —(9(1), these black holes must necessarily have 
scalar hair. In Appendix A we numerically integrate the field equations to verify whether 
c*f boundary conditions admit static, spherically symmetric black hole solutions with scalar 
hair outside the horizon. We find a one-parameter set of spherical hairy black holes, which 
can be characterized by their conserved mass. However, it turns out that for / > the 
hairy black holes are always more massive than a vacuum Schwarzschild-AdS black hole of 
the same size. Hence the singularity that develops from the spherical negative mass initial 
data defined by the instanton cannot be hidden behind an event horizon. It is therefore 
plausible that it extends all the way to the boundary, cutting off all space. 

2.4 Ten-dimensional Viewpoint 

D — 5, Af — 8 gauged supergravity is believed to be a consistent truncation of ten- 
dimensional IIB supergravity on S 5 . This means that it should be possible to lift our five- 
dimensional solution to ten dimensions. At the linearized level, the D=5 scalar fields which 
saturate the BF bound correspond to quadrupole modes on S 5 . Since in our example, the 
D=5 scalar field we retain diverges at the classical singularity, one expects that the sphere 
will become highly squashed. 

Even though it is not known how to lift a general solution of D — 5, M = 8 supergravity 
to ten dimensions, solutions that only involve the metric and scalars saturating the BF 
bound can be lifted to ten dimensions [43] . The ten-dimensional solution involves only the 
metric and the self-dual five-form. To describe them, we first introduce coordinates on S* 5 
so that the metric on the unit sphere takes the form (0 < £ < n) 

dtt 5 = d£ 2 + sin 2 £rffi 4 . (2.14) 

Letting / = e 7 ^ and A 2 = / sin 2 £ + f~ 5 cos 2 £, the full ten-dimensional metric is 

ds\ Q = Ads 2 5 + f 4 Ad£ 2 + (/A)" 1 sin 2 £dtt 4 , (2.15) 

which preserves an SO (5) symmetry of the five-sphere, as expected. The five-form is given 
by 

G 5 = Ue 5 + 6 sin £ cos Zf' 1 * df A df, (2. 16) 

where e 5 and * are the volume-form and dual in the five-dimensional solution and 

U = -3f sin 2 £ + J" 10 cos 2 £ - J" 4 - 4/~ 4 cos 2 £. (2.17) 
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In the homogeneous region of the asymptotic AdS& space, the metric can be written in 
Robertson- Walker form (12.131) . Near the singularity both the potential and the curvature 
are unimportant in the Friedmann-Lemaitre equations, and we have a(t) oc (t s — t) 1 / 4 and 
<p(t) = -(y/3/2) ln(t s - t). Therefore, over most of the S 5 near the singularity, the metric 
approaches 

ds* =sin£[-(t.-t)-V^d^ + (t 8 -^^ 

(2.18) 

Introducing a new time coordinate T — (t s — t) M / 4v/l ~°, where \x = 4a/T0 — 1, this takes a 
simple Kasner-like form 

ds 2 w = sin£[- (^j^J 2 dT2 + T^-^dai + T~ 18 ^d^ 2 + T 6 ^dQ A }. (2.19) 

The anti de Sitter space and four of the dimensions of the S 5 shrink to zero, with Kasner 
exponents of approximately 0.19 and 0.26 respectively, while the fifth dimension of the S 5 , 
labelled by £, blows up with a Kasner exponent of ~ —0.77 so the S 5 becomes spindle- 
shaped. Following the analysis of [44], we can T-dualize the £ dimension to obtain a 
homogeneous spacetime solution of type IIA string theory. The duality-invariant dilaton is 
20 — £\ Aj, where the scale factor of the i'th spatial dimension is = e A,; . Under T-duality, 
A^ — > — A^ and the string- frame Kasner exponent for the £ dimension becomes +0.77. The 
dilaton was static in the original IIB frame, but in the T-dual theory we have ~ T 77 so 
the string coupling tends to zero at the singularity. The new metric is now contracting in 
all directions, and hence qualitatively similar to the isotropic cosmological solution to the 
low energy effective action for string theory in string frame, with oc Ts for i — 1, . . . 9, 
and oc T. This background solution in type IIA theory corresponds to the solution 
to 11-dimensional M-theory in which the M-theory dimension and time form compactified 
Milne spacetime, with the other dimensions are static [4,20]. It will be very interesting to 
see whether we can find an unstable mode within the generalized AdS 5 x S 5 setup which, 
near the singularity, corresponds precisely to the collapse of the M-theory dimension and 
which could, once the appropriate boundary deformation is identified, be used to model an 
end-of-the- world brane collision in 11 dimensions, in the heterotic model. 

More generally, it is clear that the cosmological solution we focus on here is only one 
of many possible cosmologies allowed by generalized boundary conditions on AdS 5 x S 5 . 
The AdS/CFT correspondence can thus be used as a "laboratory" for the study of the 
nonperturbative counterparts (in both a' and g s ) of a large class of cosmological solutions 
to the low-energy effective actions for string theories. Every nontrivial solution possesses 
a spacelike singularity but, by mapping the theories in each case into an unstable dual 
quantum field theory, it may be possible to resolve many of these singularities and to 
describe the passage of model universes through them. 
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3 The Boundary Theory: a Double Trace Deformation of J\f = 4 
Super- Yang-Mills Theory 

In the previous section, we have seen that our bulk theory with boundary conditions (I2.10p 
allows smooth, asymptotically AdS initial data to evolve in finite time into a big crunch 
singularity that extends all the way to the boundary. Now we discuss the dual CFT coun- 
terpart of this phenomenon. First, we review that the boundary conditions f l2.10j) with 
/ > correspond to adding an unstable potential to the boundary field theory. Then we 
argue that, in this particular model, the quantum effective potential shares the property of 
being unbounded below. 

For the usual a = boundary conditions on the bulk scalars, the dual field theory is 
M = 4 super Yang-Mills theory. The bulk scalars that saturate the BF bound in AdS 
correspond in the gauge theory to the operators cTr[$ l $ J ' — (l/6)5 y $ 2 ], where are the 
six scalars in M = 4 super Yang-Mills and c is a normalization factor that will be fixed 
momentarily. The S'0(5)-invariant bulk scalar cp that we have kept in (12.61) couples to the 
operator 



O = cTr 



(3.1) 



According to the AdS/CFT correspondence, this means the following. In the asymptotic 
behavior (I2.8p . the function a of the coordinates along the boundary plays the role of a 
source for O in the field theory: the field theory action has a term J d 4 x a{x)0(x). The 
usual boundary conditions set this source to zero, meaning that the boundary field theory 
is the undeformed M = 4 super Yang- Mills theory. On the other hand, the function f3 plays 
the role of the expectation value of O in the field theory; different (3 correspond to different 
quantum states of the boundary theoryjfl 

In general, imposing nontrivial boundary conditions a{j3) in the bulk corresponds to 
adding a multi-trace interaction W(0) to the CFT action, such that after formally replacing 
O by its expectation value j3 one has [13,47] 

SW . , 

a = - w (3.2) 

Adding a source term to the action can be considered as a special case, where W is a 
single-trace interaction linear in O. The boundary conditions (12.101) that we have adopted 
correspond to adding a double trace term to the field theory action 

S = S - W{0) = Sv + f -J(D 2 . (3.3) 

The operator O has dimension two, so the extra term is marginal and preserves conformal 
invariance, at least classically (we shall see that conformal invariance is broken quantum 
mechanically) . 



8 For a detailed discussion of the relation between the bulk field and the Yang-Mills operator, see for 
instance [45]. For issues specific to AdS/CFT in Lorentzian signature, see for instance [46]. 
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In the previous section we have taken the constant / to be small and positive in the 
bulk. The term we have added to the CFT action, therefore, corresponds to a negative 
potential. Since the energy associated with the asymptotic time translation in the bulk 
can be negative, the dual field theory should also admit negative energy states and have 
a spectrum unbounded below. This shows that the usual vacuum must be unstable, and 
that there are (nongravitational) instantons which describe its decay. After the tunneling, 
the field rolls down the potential and becomes infinite in finite time. This provides a 
qualitative dual explanation for the fact that the function (3 of the asymptotic bulk solution 
( 12. 8p diverges as t — > tt/2, when the big crunch singularity hits the boundary. Since f3 is 
interpreted as the expectation value of O in the dual CFT, this shows that to leading order 
in 1/N, (O) diverges in finite timej^l 

So the big crunch spacetime in the bulk theory corresponds in the boundary theory 
to an operator rolling down an unbounded potential in finite time. It is important to 
know whether quantum corrections preserve the unbounded nature of the potential in the 
boundary theory. While this was unclear for the AdS4 model that was the main focus 
of earlier work [10,11,48,49], we shall now argue that for our model we indeed have an 
unbounded quantum effective potential. 

For this purpose, we first briefly summarize the renormalization properties of a double 
trace deformation (13. 3p of M = 4 super- Yang-Mills theory; a more detailed discussion can 
be found in Appendix B. As explained in [13], the computation of amplitudes at order f 2 
involves matrix elements of 

£ J d 4 xd 4 y0 2 (x)0 2 (y). (3.4) 

From conformal invariance, (O(x)O(y)} = v/\x — y\ 4 on flat H 4 (where the constant v 
depends on the normalization factor c in (13.1 ft as well as on N). This leads to a short 
distance divergence that renormalizes / and survives in the large N limit: 

^ J d 4 xd A yO(x)0(y)(0(x)0(y)) ~ n 2 f 2 vlnA j d 4 x0 2 (x), (3.5) 

with A an ultraviolet cutoff. This leads to a one-loop beta function for /, which does 
not receive higher loop corrections in the large iV limit [13]. As in [13], we now fix the 
normalization constant c (and thus v) by demanding that the beta function coefficient 
should be one. At least for small 't Hooft coupling, it is easy to see that v ~ c 2 N 2 for 
large N\ therefore c = a/N with a a numerical constant. The coupling / can then be kept 
fixed in the large A" limit. The existence of a non- vanishing beta function means that the 
conformal invariance of Af = 4 super- Yang-Mills theory is broken quantum mechanically 
by the double trace deformation. 

9 In fact, we shall see later that when we consider our theory on R x S 3 , the expectation value of O 
in a state described by a wavepacket diverges well before the center of the wavepacket reaches infinity, so 
in our deformed theory on a finite volume space it is inappropriate to phrase the dynamics in terms of 
expectation values. What happens in that case is that the bulk of the wavepacket reaches infinity in finite 
time. 
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We will be interested in an approximation to the quantum effective action that is valid 
for a large range of field values, in particular for large field values. An appropriate frame- 
work is that of [50] , where the standard Feynman diagram expansion is resummed and the 
theory is organized in a derivative expansion and an expansion in the number of loops (see 
Appendix B). The one-loop effective potential is given by [51] 

f u 9 f 2 2 \n(OI^ 2 ) 
V(0) = --^0 2 + ^ V 7 , (3.6) 

where \x is a renormalization scale, f ^ a renormalized coupling, and the counterterms have 
been chosen such that there is no constant and no linear term in O, and such that 

W) = -^/A (3.7) 

The renormalization group equation can be obtained by demanding that V(0) be indepen- 
dent of \i: 

^ = ~fl (3-8) 

which shows that the normalization of O implicit in (I3.6P is indeed such that the beta 
function coefficient is one. In (13. 8p . we have ignored a contribution from d/dfj, hitting the 
f 2 in the second term of (13.61) . which is justified as long as \f^ ln(C/yU 2 )| <C 1. Equation 
(13. 8p is solved by 

fa = (3.9) 

ln(///M)' 1 ' 

with M an arbitrary scale (this implements dimensional transmutation). Choosing /i 2 = O, 
i.e., the renormalization scale is set by the value of the field O, the Coleman- Weinberg 
potential can then be written as 

O 2 

V(0) = (3.10) 

V ; ln(0/M 2 ) V ; 

Now suppose that for some value C , the coupling is small, 

< f O0 « 1, (3.11) 

then M 2 < Oq and (13.101) is trustworthy (i.e., higher order corrections can be ignored) for 
any O such that O > Oq. As a result, we can conclude that in this case 

V(0) -> -oo for O -> oo. (3.12) 

Other quantum corrections are small in the large O regime we will be interested in, as 
described at the end of Appendix B. 

As was shown in [13], these renormalization properties have a precise counterpart in the 
bulk theory. First, define a new coordinate 0<n<lbyr = 2u/ (1 — u 2 ). In terms of this 
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coordinate, the conformal boundary of AdS is at u — 1. Near a point on the boundary, the 
metric takes the form 

ds 2 = ^^(-dt 2 + dz 2 + dx* dx*), (3.13) 
z l 

where z — 1 — w ~ 1/r and x l replace the coordinates of the three-sphere (which is approx- 
imately flat when zooming in on one point) [52]. With the boundary condition a = ff3, the 
behavior of the field near the boundary z = is 

(f) = p z 2 (-f\nz+ 1). (3.14) 

It was shown in [53] that the near-boundary (small z) region of the bulk theory corresponds 
to the UV of the dual field theory. Therefore, to study this UV regime, we introduce a new 
coordinate 

~z = Z - (3.15) 

with e < 1; this new coordinate is well-suited to studying the near-boundary (small z) 
region of AdS. In terms of z, the boundary behavior (13.141) reads 

= f3z 2 (-f\nz + l) (3.16) 

with 

/= t4i^ (3 ' 17) 

and $ = (1 — f\ne)/3. Interpreting e as a ratio of renormalization scales, e = fi/fi, and 
interpreting / = and / = fp, as the coupling defined at the scales \i and fl, respectively, 
(13.171) implies the following relation between the couplings at different scales: 

U = i FT7 /~n> (3.18) 

1 - ln(/x//x) 

which is consistent with the renormalization group equation (13.81) . The reason that this bulk 
computation (valid for large 't Hooft coupling) agrees with the perturbative field theory 
computation (valid for small 't Hooft coupling) is the fact that the beta function is one-loop 
exact for large N [13] 

In what follows we will mostly concentrate on the steepest negative direction of the 
effective potential. Fluctuations in orthogonal directions in field space acquire a positive 



10 To relate our version of the above argument to Witten's, given in Ref. [13], note that our conventions 
are related to his by z = rw,a = —aw,f — —fw- In the original argument, the boundary condition is 
written as 

cj) = /3z 2 (-/ln(^z) + 1), (3.19) 

where /i is an arbitrary scale introduced to define the logarithm. One can choose a different mass scale fl 
if one "renormalizes" the field (3 and the coupling / in such a way that the bulk field <fi is left invariant: 

P{-fMpz) + 1] = /?[-/ ln(,iz) + 1], (3.20) 

which implies the relation (|3.18|) between the couplings at different scales. 
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mass and we will see these are suppressed. For the SO(5)-invariant operator we consider, 
the most unstable direction comes from the — $f term in (13.101) (see (13.11) ). We focus on 
the dynamics of $1 as it rolls along a fixed direction in su(N): 

$i(x) = (j){x)U (3.21) 

with U a constant Hermitian matrix satisfying Trf/ 2 = 1, so that is a canonically nor- 
malized scalar field. The Coleman- Weinberg potential (13.101) for this scalar is then given 
by 

VW = = "T ^ ~~i • (3 - 22) 

where 

2a 2 9 M 2 , A , 

with a the numerical constant implicitly defined after (13.51) . 



4 Unbounded Potentials, Self- Adjoint Extensions and Ultra-Locality 

We have seen that our field theory description involves the potential (I3.22j) . which is un- 
bounded below. This implies that the quantum field theory has no ground state. While 
such unstable theories are usually considered unphysical, we want to explore whether quan- 
tum mechanical evolution can be defined for them in a consistent way. A clue is provided 
by quantum mechanics (as opposed to quantum field theory) with unbounded potentials. 
As we shall review momentarily, if a potential allows a wavepacket to move off to infinity in 
finite time, one can nevertheless define unitary quantum evolution by imposing appropriate 
boundary conditions at infinity. Technically, one restricts the domain of allowed wave- 
functions to those on which the Hamiltonian is self-adjoint - this is called a "self-adjoint 
extension" . In this section, we begin to develop the case that the quantum field theory we 
are interested in should also possess a self-adjoint extension. In Subsection 14.11 we review 
unbounded potentials in quantum mechanics. In Subsection 14.21 we show explicitly how the 
field theory dynamics become "ultralocal" as the singularity is approached, thus making it 
plausible that the quantum field is described by an independent set of identical quantum 
mechanical systems, one for each spatial point. 



4.1 Quantum Mechanics in a — Ax p /4 Potential for p > 2 

We are interested in a theory of a scalar field which is classically conformal invariant but 
unstable. Let us first emphasize the generality of this setup, in a holographic context. We 
consider a classically conformal-invariant scalar field theory on a d— dimensional conformal 
boundary of the form 1R x S^ 1 where 1R is time. The action is 

J d d x^f d (~(<90) 2 - £A0 2 + i V p ) , (4-1) 
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where Rd is the Ricci scalar, and the coupling A p is an arbitrary constant. Conformal 
invariance requires p = 2 + 4/(g? — 2) and ^ = (d — 2)/ (4(d — 1)). For any d > 2, p is 
greater than 2 and provided X p is positive, will run to infinity in a finite time. When 
focusing on the behavior near the singularity, where the <fr p term dominates in the potential, 
we shall neglect the Rd4> 2 term. 

The simple fact that the spatial volume of the conformal boundary is finite is very impor- 
tant for the quantum behavior of the boundary theory. Consider the quantum description 

of the homogeneous mode (f> of the scalar field. Its kinetic term in the action is Vd-\ J dt\^> 
where the volume of space, Va-i, acts as the "mass" of <fi. In the infinite volume limit, this 
"mass" becomes infinite, and <fi undergoes no quantum spreading: it becomes a classical 
variable. (This is the essential reason why spontaneous symmetry breaking in possible in 
quantum field theory but impossible in quantum mechanics). When V^-i is finite, as here, 
the homogeneous mode (f) undergoes quantum spreading. It is convenient to canonically 
normalize the homogeneous mode, setting x = (Vd-i)^(po and A = (V^_i) 1 ~ 3 '/ 2 Ap, which is 
constant (for now, we are ignoring the running of the coupling constant). We then have a 
unit mass quantum mechanical particle with coordinate x and potential 

V(x) = ~\x p . (4.2) 

In this subsection, we summarize the operator approach to the quantum mechanics of such 
potentials [14] as reviewed in [15] In Subsection 15.21 we shall implement the self-adjoint 
extension using complex solutions to the classical equations of motion, which will be useful 
for describing the evolution of Gaussian wavepackets. 

A classical particle rolling down the potential (14.21) . with p > 2, reaches infinity in 
finite time. The same is true for a quantum mechanical wavepacket, and this at first 
sight appears to lead to a loss of probability, i.e., to non-unitary evolution. However, if 
a self-adjoint Hamiltonian could be defined for this system, unitary quantum mechanical 
evolution would be guaranteed. As we shall now review, this can be done by carefully 
specifying an appropriate domain for the Hamiltonian 

/t = - l -*3 + VM. (4.3) 

Since the WKB approximation becomes increasingly accurate at large x, we can use it to 
study the generic behavior of energy eigenfunctions there. The two WKB wavefunctions 
for fixed energy E are proportional to 

X %{x) = [2(E + AX74)]- 1 / 4 exp (±i j ^2(E + \yP/4)dy) , (4.4) 

where the lower limit of the integral may be chosen arbitrarily. 



11 For a related recent discussion, see [54]. 
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The Hamiltonian is self-adjoint if for any wavefunctions 4>i, 4>2 in its domain, (H<f>i, <pz) — 
(<fii, Hfa)- Using integration by parts, one sees that this is equivalent to 



dx dx 



(4.5) 



This can be arranged if for each energy E we select the linear combination of the two WKB 
wavefunctions which behaves like 



\/2\r p / 2+1 \ 

ip%(x) ~ x~ p/4 cos — — + a (4.6) 

y p+2 J 

at large x, where a is an arbitrary constant phase. (The angle < a < tt labels a one- 
parameter family of inequivalent self-adjoint Hamiltonians.) If E is positive, for example, 
we set 

r E (x) = [2(E + Aa;V4)]- 1 /4 C0S Qf ^ 2 (E + \yP/4)dy + <p%j (4.7) 



with 







y/2(E + Ayf/4) - <fy, (4. 



which tends to the required form at large x. A similar construction can be given for negative 
E [15]. 

As a consequence of the fact that every energy eigenfunction tends to the same, energy- 
independent form (I4.6p . equation (14. 5 p is satisfied if 0i = ip^ and 2 = ^e 1 with the same 
value of a. The domain of the "self-adjoint extension" of H labelled by a can now be 
defined as all wavefunctions 

(f>(x) = J dE4>(E)^ E (x) (4.9) 

that satisfy J dE\cj)(E)\ 2 < oo, so is square integrable, and J dEE 2 \(p(E)\ 2 < oo, so 
that H(p is also square integrable, i.e., it is also a normalizable wavefunction. Under these 
conditions, the inner product (0i, Hfa) makes sense and, as we have just checked, it equals 
(if0i,0 2 ) so that H is self-adjoint. 

One can interpret the parameter a as follows [15]. If one placed a "brick wall" at large 
x, it would force the energy eigenfunctions to vanish there. However, since for finite E the 
de Broglie wavelength becomes independent of energy at large x, displacing the wall by a 
half-integral number of de Broglie wavelengths would, in this regime, have no effect. Hence 
only the location of the brick wall modulo an integer number of half-wavelengths matters 
physically, and this is the information contained in the phase a. 

Physically, one can interpret these self-adjoint extensions as follows. A right-moving 
wavepacket that moves to infinity is always accompanied by a left-moving "reflected" 
wavepacket that runs back up the hill. The time it takes for a right-moving wavepacket 
to run to infinity and for the left-moving wavepacket to run back up the potential hill can 
be shown to be the same as a classical particle would take to fall to infinity and climb up 
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the potential again after being reflected at infinity. If the potential is bounded for x < 0, 
the Hamiltonian allows a continuum of scattering, states. It also has an infinite number of 
bound states with quantized (negative) energies] 12 ! the values of which depend on a. How- 
ever, for the scattering problem - the bounce off the singularity - which we study, to a very 
good approximation the phase a enters only as an overall phase in the final wavefunction 
and hence has no physical consequence. 

In Section 6, we shall compute the Schrodinger wavefunctional ^ (0(x)) for the full 
quantum field, decomposed into its homogeneous and inhomogeneous parts, 0(x) = + 
<50(x). The homogeneous part 0, discussed in Section 5, behaves like the coordinate x 
considered in this section. In Subsection 15.21 we take the initial wavefunction for to be 
a Gaussian with a certain spread, specified in Subsection 15.41 The inhomogeneous part 
<50(x) is well- described at early times in terms of its Fourier modes, describing a set of 
harmonic oscillators which we shall take to be in their incoming ground state. We shall 
argue in Subsection 6.6 that for all but a narrow band of final values of centered on the 
real classical solution, we can treat the inhomogeneous modes 50(x) to quadratic order 
while ignoring their backreaction on 0. Therefore, we conclude, the final wavefunction may 
be calculated reliably, using perturbation theory in <50, for a range of values of away 
from the real classical solution but carrying most of the probability. Provided we are in a 
certain parameter range, to be detailed in Subsection 6.6, the system undergoes a bounce 
and returns close to its initial configuration, with little backreaction from particle creation. 

4.2 Ultralocality and Self-Adjoint Extensions in Quantum Field Theory 

To address the question whether these methods can be extended to quantum field theory, 
we shall first consider a simplified model that shares the same finite time singularity, namely 
the quantum field theory of a single scalar field with a negative quartic potential. In the 
context of the actual dual quantum field theory, this approximation amounts to concentrat- 
ing on the scalar that parameterizes the steepest negative direction: the gauge-invariant 
magnitude of $i. We shall argue that the method of self-adjoint extensions indeed applies 
to this quantum field theory. The main argument is that near the singularity, the scalar 
field evolution becomes ultra-local: spatial gradients become unimportant and the quantum 
field theory can be thought of as a collection of independent, identical quantum mechanical 
systems, one for each point in space. As a consequence, the self-adjoint extension method 
can be implemented point by point. We assume that at every spatial point, the same value 
of a enters the self-adjoint boundary condition: this is necessary if the theory is to respect 
the symmetries of the S 3 . 

It is worth elaborating on this point. A priori, one could imagine self-adjoint extension 
phases ct(x) that depend on the position x on the three-sphere. Consider the 50(4) rotation 

12 We note in passing that there is another approach to quantum mechanics in potentials which are 
unbounded below, based on a PT symmetry [55]. This approach is motivated by analytic continuation 
from the harmonic oscillator and it results in a positive energy spectrum. In contrast, for all self-adjoint 
extensions described above, there is an infinite set of negative energies. Since we know the bulk theory has 
negative energy solutions (the instanton discussed in Section 2 provides one example), the approach using 
self-adjoint extensions seems more appropriate. 
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group of the S 3 . When choosing a theory, for instance by specifying a Lagrangian, it is 
considered natural to preserve this symmetry group; for instance, the Lagrangian density we 
are about to consider does not have explicit x-dependence. On the other hand, choosing a 
symmetric state is often considered an undesirable fine-tuning. Our point is that the choice 
of self-adjoint extension is part of choosing a theory, not a state. Indeed, a Hamiltonian 
with a steep unbounded potential is not defined as an operator until its domain is specified. 
More physically, the corresponding theory is not defined until its space of states is specified. 
Therefore, we consider it very natural to choose the self-adjoint extension phase a to be 
independent of x, which is the choice we shall make in the present paper. Having said this, 
it would be perfectly legitimate (albeit less natural) to consider other choices as well (just 
as one could consider Lagrangian densities with explicit x-dependence). 

To begin with, we shall study the classical, conformally- invariant field theory onlRxS' 3 , 

C = - 1 -(d4>r + ^-^R^, (4.10) 

where -R4 is the Ricci scalar. We are interested in studying generic solutions, i.e., those 
specified by the appropriate number of with arbitrary functions of space, in the vicinity 
of the singularity, i.e^ the locus where (f) reaches infinity, which we shall assume to be a 
space-like surface £0 The field equation is 

□0 = -A 4 3 + -R A( j). (4.11) 
o 

The curvature of the three-sphere plays an important role in determining the initial con- 
ditions for the scalar field, but as the field runs down towards the singularity the spatial 
curvature rapidly becomes irrelevant. Since we are concerned with the behavior of the field 
near the singularity, in this section we shall ignore the last term in (14.111) . 

It is also convenient to change variables to x — (2/A4)2 _1 so that the equation of 
motion assumes the form 

X d 2 X - 2(d X ) 2 = 2, (4.12) 

and the zero-energy, homogeneous background solution is just x — t* ~ ^ f° r t < t*, with 
t* an arbitrary constant. 

We would like to construct the generic, spatially inhomogeneous solution to (14.121) as a 
Taylor series in the proper time away from the singularity, which, as mentioned above, we 

13 The form of S will obviously depend on the specific solution one considers, i.e., on the initial conditions 
for the field <f>. One may wonder whether £ will be a smooth, spacelike surface for reasonable initial 
conditions, and in particular for those assumed in the remainder of the paper, where the inhomogeneous 
field modes start out in the adiabatic vacuum. These quantum mechanical initial conditions generate 
a classical solution by the phenomenon of "mode freezing," essentially in the same way that classical 
perturbations are generated from quantum fluctuations in inflationary cosmology. For those solutions, we 
shall verify in Subsection 6.6 that the surface £ is indeed smooth and space-like, justifying the assumption 
we are presently making. Note that this conclusion is not in contradiction with the conclusion we shall 
reach in the present section, namely that the field will evolve in an ultra-local way near the singularity. 
While ultralocality implies that field fluctuations 8<j) grow as the singularity is approached, it does not cause 
wild fluctuations in the surface S, as we shall show in Subsection 6.6. 
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Figure 5: The singular hypersurface E, assumed spacelike, upon which <\> is infinite and \ 
is zero. 



assume to be a spacelike surface E. The coordinates of the embedding Minkowski spacetime 
(recall, we are neglecting the space curvature of the S 3 ) are t, x; the locus of the singular 
hypersurface E is t = t s (x s ), x = x s . As long as E is not too strongly curved, for every 
spacetime point x M there is a unique corresponding point x^ on E, connected to x M by a 
geodesic normal to E. The unit normal to E is 

n» = (7, 7Vt s ) , 7 = . (4.13) 

where V indicates the gradient with respect to x s . It is then natural to change coordinates 
from (t, x) to r, the proper time from E, and x s , the corresponding point on E, given by 
x M — x^ = n^r, or 

t = i s (x s )+ 7 r, x = x s + 7 rVt s , (4.14) 

where r = yjt — t s ) 2 — (x — x s ) 2 is (the absolute value of) the proper time from E. 

To construct the scalar field equation we need the metric in the new coordinates. From 
the Minkowski line element in t, x, and (14.141) . we find the line element in (r, x s ) coordinates 
is —dr 2 + hijdx l s dxi, where 

= +2K ij r + K ik h% ) K l jT 2 
h if = $ij ~ dit s djt s , K i:j = ^didjt s . (4.15) 

Here, t s (x s ) is the locus of the singular hypersurface, whose metric is and inverse metric 
is /i/Q\ . All space derivatives di are taken with respect to x s , at fixed r. The matrix Kij is 
the extrinsic curvature (or second fundamental form) of the singular hypersurface E. (Since 
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the embedding space is flat, the intrinsic curvature of £ can also be expressed in terms of 
Kij, by the Gauss relation Rj kl (E) = Kj k K\ — KjiKl, and is thus not an independent 
quantity.) 

In these coordinates, the field equation (I4.12p reads 

2x 2 T - xh'HhhAr ~ 2h ij d iX d jX + ^"s^U^,^ = 2, (4.16) 

where h is the determinant and h l i the inverse of hy. By substituting (14.151) into (14.161) . 
we obtain an asymptotic expansion for x, 

1 1 1 1 3 7 1 

a 6 lg v i 2j 4 V 3 18 54 1 6 s ; 



+lnr 



5 V 3 3 z 9 1 27 1 11 

r 5 + J D 6 r e + ..., (4.17) 
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where K\ = h^Kij = h 2d T hi\ T= o, K 2 = Tr(K 2 ), and so on, and traces are taken using 
the metric hV Q y The higher coefficients C n (x s ), n > 6, are all determined in terms of t s (x s ), 
and the higher coeffients D n (x. s ), n > 6, are determined in terms of the second arbitrary 
function p(x s ). The reason for the entrance of the logarithm, and the second arbitrary 
function, at order r 5 is seen by examining the first two terms of (14.161) . at order r 4 . At this 
order, it is necessary for consistency to introduce a term proportional to r 5 In r into the 
series expansion, and the coefficient of r 5 is then left undetermined. 

Let us comment on the general form of the full nonlinear solution (14.171) . First, if we 
re-express r in terms of the global time coordinate, 

r= |t-t s (x s )|+0((Vt s ) 2 ), (4.18) 

(14.171) can be viewed as an expansion in powers of Vt s as well as in (t — t s ) V s . However, if 
we linearize in the two arbitrary functions t s (x) and p(x), we obtain the solution 

X (t, x) = -t + t s (x) + Vv 2 t,(x) - ^-t 4 V%(x) + . . . - 1 5 + . . . , (4.19) 

6 24 10 

which is actually completely regular at t = 0. Within this linearized approximation, the two 
arbitrary functions are easily interpreted. We decompose the scalar field into a homogeneous 
component <ft, treated nonlinearly, and an inhomogeneous component 50(x) which we treat 
only at the linearized level. Then £ s (x) is the time delay to or from the singularity, given in 

linear theory by the limit of 8(j)/<f> as t tends to zero. And p(x) is the linearized perturbation 
in the Hamiltonian density, evaluated at the singularity: setting 

6H = %5<P + VS)W = dt[ 5 4)f, (4-20) 
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we find that as t — > 0, 57i tends to p(x). Note that spatial gradients do not appear 
in SH at linearized order because the background solution <fi = — ^/A^t -1 is spatially 
homogeneous. 

Second, (I4.17P is clearly an expansion in rV s , so that gradients become less and less 
significant in the dynamical evolution as the singularity r = is approached: the evolution 
of x becomes ultralocal in this limit. To see this, consider the "velocity-dominated" version 
of the equation of motion (I4.12p . i.e., the equation of motion with the gradient terms 
omitted. The solution of the velocity-dominated equations of motion is 

X (0) =t s (x)-t + ^(t s (x)-t) 5 . (4.21) 

Therefore, we see that the full solutions (14.171) differ from the velocity-dominated ones only 
by terms that are unimportant near the singularity, so that the gradient terms are indeed 
dynamically unimportant there. 

Finally, note that the non-analytic terms in (I4.17p . involving lnr, enter only at quadratic 
and higher order in the perturbations. If these nonlinear terms are constructed perturba- 
tively by expanding in the linearized perturbation, then no new functions enter: the nonlin- 
ear terms are completely dictated by the equations of motion. Hence, a matching condition 
within linear theory completely determines the matching of the full nonlinear solution, as 
long as the latter is well- described by linear perturbation theory. 

To summarise, we have exhibited a generic class of classical solutions which are singular 
on an arbitrary spacelike hypersurface S. The solutions are constructed as an expansion 
in proper time r times spatial gradients: if we follow a solution with spatial structure on 
some wavenumber k, gradients become unimportant in the evolution when \kr\ falls below 
unity. At this point, the evolution becomes ultralocal so that the field values at different 
spatial points decouple. This behavior is relevant to the full quantum theory, because the 
semiclassical expansion becomes exact as the singularity is approached, so that the quantum 
wavefunction can be constructed from classical solutions. 

In contructing the quantum theory, it is then natural to use the same self-adjoint ex- 
tension (or boundary condition at large 0) at every spatial point, to evolve the quantum 
field across the singularity (see the discussion of this point at the beginning of this subsec- 
tion). We shall implement the self-adjoint extension using the method of images, so the full 
wavefunction will be constructed from two classical solutions with one, roughly speaking, 
starting out "behind" the singularity. A key question is whether the UV regime of the 
quantum field theory can be consistently treated, or whether the presence of structure on 
smaller and smaller scales will lead to divergences. The remainder of this paper will largely 
be devoted to tackling this question. 

In fact, we shall not need the full nonlinear solutions given above. For the most part, 
our calculations will be done in linearized perturbation theory around the homogeneous 
background. We start from a state in which the homogeneous mode is described by a 
localized wavepacket, up the hill, and the inhomogeneous modes are in their adiabatic 
vacuum. We shall show in Subsection 6.6 that effects nonlinear in the inhomogeneities are 
consistently small (in a certain parameter range, to be detailed there). The ingoing and 
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outgoing states are most conveniently described using Fourier modes for 50 and, since we 
work only at the linearized level, this description is valid throughout. The k 2 term in the 
equation of motion for these terms becomes insignificant near the singularity, reflecting the 
phenomenon of ultralocality possessed by the full theory. 

5 Quantum Evolution of the Homogeneous Component 

We want to compute the Schrodinger wavefunctional \l/ (t, 0(x)) for the full quantum field 
0(x) in the semiclassical expansion. We write the field as the sum of its homogeneous 
and inhomogeneous parts, = + 50(x), and assume that each Fourier mode of the 
inhomogeneous component starts out in its quantum ground state. As we shall see later, 
the inhomogeneous quantum fluctuations are suppressed by powers of A^), so it makes sense 
to first calculate the wavefunction for homogeneous component 0, with the inhomogeneous 
modes ignored. 

As emphasized in Subsection I4.1[ is quantum mechanical because the boundary the- 
ory lives on a finite space. This simple fact plays an essential role: the quantum spread 
in means that the bulk of the wavefunction stays away from large 0. In the semiclas- 
sical expansion, this will translate into the fact that the relevant classical trajectories are 
generically complex and generically avoid the singularity at infinity in the complex 0-plane. 

In this section, we show how the self-adjoint extension described in Subsection 14.11 can 
be implemented in the semiclassical expansion. Since the leading term in this expansion 
becomes exact near the singularity, the calculation should accurately capture the physics 
there. We impose the required boundary condition on the wavefunction at infinite by the 
method of images, using a reflection symmetry of the Hamiltonian to satisfy the boundary 
condition with the sum of two semiclassical wavef unctions. The resulting wavefunction 
bounces from infinity and returns close to its starting point. The quantum spread in the 
wavefunction is represented by a small imaginary part in each of the relevant complex 
classical solutions which is, for generic final values of 0, nonzero at all times and which 
dominates at times corresponding to the classical encounter with the singularity. 

Having determined the wavefunction \l/(0) in the approximation where the inhomoge- 
neous modes were ignored, the next step is to construct the full wavefunction \l/ (tf, 0, 50(x)) 
at the final times of interest, ignoring the backreaction of 50 on 0. We again employ the 
semiclassical expansion, finding the relevant (spatially inhomogeneous) complex classical 
solution which incorporates the initial condition that the inhomogeneous fluctuation modes 
of 50 are in their incoming adiabatic ground state. Working to quadratic order in 50 in 
the action and ignoring backreaction, the solution for 50 is determined from the relevant 
linearized field equation. The small imaginary part in the complex solution for the back- 
ground 0, representing its quantum spread, now plays an essential role: it provides a UV 
cutoff on the production of high k modes. Including this effect, we find that, within a cer- 
tain parameter regime, for generic final values for particle production is parametrically 
suppressed and backreaction is under control. 

The combined inhomogeneous classical solution yields \l/ (tf, 0, 50(x)) , in a consistent 
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approximation scheme, for values of carrying most of the final probability. The approx- 
imation breaks down for values of very close to the real classical solution corresponding 
to the initial center and momentum of the Gaussian wavepacket. As the final value for 
is taken near to this classically predicted value, the imaginary part of the relevant complex 
classical solution disappears. If backreaction is neglected, the quantum particle production 
of inhomogeneous modes diverges in this limit, indicating that our approximation breaks 
down. However, the relevant band of final yielding strong backreaction is parametrically 
small within our chosen parameter regime. For all other values of 0, our calculation of the 
joint probability density (tf,(f), 50(x) | 2 appears to be reliable. Since we are working in 
a theory where unitarity has been built in, we conclude that the most probable outcome is 
a bounce, with parametrically small particle production. 

We begin this section with a review of the use of complex classical solutions to construct 
the solution of the time-dependent Schrodinger equation in a semiclassical expansion. 



5.1 Complex Classical Solutions and Quantum Mechanics 

Our discussion of the evolution of the quantum field will rest upon the use of complex 
classical solutions. These are used to obtain the first semiclassical approximation to the 
Schrodinger wavefunction, which then forms the basis of the semiclassical expansion. The 
method is ideally suited to the present problem because the first semiclassical approximation 
becomes exact as the singularity is approached. 

Consider a quantum mechanical particle with Hamiltonian H(x,p) = \p 2 + V(x). At 
the initial time t = tj we prepare the system in a Gaussian wavepacket centred on x = x c 
and with momentum p = p c . Our initial condition may be written as 

x .pL x c .p c L 

Vi — = V i = a, t = ti, (5.1) 

because this equation, when regarded as an eigenvalue equation with p = —ihd/dx, yields 
for the wavefunction 

y{x,U)~ e l —e (5.2) 

The corresponding probability density is a Gaussian with variance Ax 2 = L 2 . Equation 
(15. ip can also be written as 

x . p x c . p c , 

"7 1" l ~r- — "7 1" l ic- ? (5.3) 

Ax Ap Ax Ap v ' 

because for a Gaussian, AxAp = h/2, saturating the uncertainty relation. 

For a given choice of A set of Gaussian wavepackets with x c and p c running 

over all real values can be thought of as providing a smeared description of classical phase 
space (x,p), with the minimal "smearing area" AxAp. 

The normalized states labeled by x c and p c , 



■p c x (x-x c ) 2 



= T^TT^-e-^i^. (5.4) 



(2ttL 2 ) 
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are sometimes called "coherent squeezed states" and the complex number a defined in (15. ip 
is the coherent state parameter. For given L, these states form an overcomplete basis: they 
are not orthogonal, but any state can be expressed as a linear combination of them. To see 
this, note that 

^*w>)*^>') = ~ (5-5) 

gives the resolution of the identity operator. An arbitrary wavefunction *f>(x) can be ex- 
pressed as a linear combination of the ^ cPc by multiplying each side of (15.51) by ty(x') and 
integrating over x'. Clearly, if we determine the time evolution of Gaussian wavepackets, 
for some L, we will have determined the time evolution of any state. 

For a harmonic oscillator with frequency to, i.e., V = \oj 2 x 2 , the semiclassical approxi- 
mation is exact. Furthermore, in this case there is a natural choice for L given by the width 
of the ground state L 2 = h/(2u). The left hand side of (15.11) is then just the annihilation 
operator a = \Ju /(2h) (x + ip/u) obeying [a, aH = 1, and in fact this was the reason for 
adopting the normalization of the coherent state parameter a in (15. 11) . The initial condition 
corresponding to the oscillator being in its ground state is simply a = 0: we shall adopt 
this initial condition for the inhomogeneous modes of the scalar field on the S 3 which, at 
times well away from the bounce, are accurately described as harmonic oscillators with 
adiabatically varying frequency u. More generally, one can consider initial and final states 
which are eigenstates of a with eigenvalue a, just as in (15.11) . These eigenstates are given 
by \a) oc e aat |0). The transition amplitude between such states (/3, out|a, in) provides the 
generating function of the S-matrix between arbitrary n-particle Fock states. [56] Therefore 
the formalism of complex classical solutions extends naturally to a computation of the full 
S-matrix, in the semiclassical expansion. 

The semiclassical expansion may be derived as a formal expansion in powers of h: sub- 
stituting 

V(x,t) = A{x,t)e iS ^ t),h (5.6) 
into the time-dependent Schrodinger equation 

( h 2 d 2 \ 

^H-w+^'r- (57) 

at leading order in h we obtain the Hamilton- Jacobi equation for S, namely 

OS ( 1 fdS\ 2 Jr/ \ TT , OS, 

This equation is solved by the classical action, i.e., S = Sci being the action evaluated on 
the appropriate classical path xci- Since our initial condition (15. ip is complex, we shall in 
general be interested in complex solutions of the classical equations of motion. The classical 
action is 

S C i(x,t) = / dt' (pciici - H(x C i,Pci)) (t') + B(xi,pi) (5.9) 
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for a classical path satisfying the initial condition (15. ip at time U and the final condition 
xci — x at the final time t. (We often abbreviate the final coordinate Xf as x and the 
final time tf as t: the context should make our meaning clear.) It is necessary to add the 
boundary term B to ensure that the stationarity of the action SS = under all variations 
5x(t) respecting the given boundary conditions is equivalent to the classical equations of 
motion. Our initial condition (15.11) implies 5xi + 2i5piL 2 /h = 0, and the correct boundary 
term is found to be 

B(x uPl ) = ^. (5.10) 

Now consider the value of the classical action Sci when the final time t and the final 
position x are varied, t —> t + 5t, x — > x + 8x. The variation in the corresponding classical 
path is Xci{t) — > Xci{t) + Sxci(t), and the variation of the classical action is 

SSci = St(pcixci - H(xci,pci))(t) + \pciSxci\l. + SB (5.11) 

where we used integration by parts and Hamilton's equations for the classical solution, 
x ci — {dH/dp)ci, Pa — ~{dH/dx)ci- The value of the final coordinate at the final time 
t + St has the first order variation 

Sx = Sxci (t + St) = Sxci (t) + xciSt, (5.12) 

hence 

SSci = -HSt + pSx - piSxi + SB = -HSt + pSx, (5.13) 

because B was chosen precisely to ensure that SB — pSxci vanishes at t — t% for our chosen 
initial condition. Thus we have 

dS C i , . dS C i f s 

- W = -H(x,p), - w =p, (5.14) 

and the wavefunction (15. 6p . with S = Sci, satisfies both the initial condition (15. ip and the 
time-dependent Schrodinger equation (15. 8p to leading order in h. 

For completeness, let us give the equation for the time evolution of the prefactor A(x, t): 

dA OA _ _ldp ihd 2 A 

~dt +P dx~~ ~2d~x + 2dx^- (5 - 15) 

With the neglect of the last, higher order term, this equation is straightforwardly solved by 
integrating A along the classical trajectories of the system. For example, if A is independent 
of t, then the solution is just A oc p~? . 

As a simple illustration, consider a free particle, with H = \p 2 ■ In this case the La- 
grangian is quadratic and the second variation of S with respect to x, i.e. (dp/dx), depends 
only on t. Here, as can be seen from (15. 15ft . it is consistent to set (dA/dx) zero. The IE- 
dependence of the wavefunction then arises solely from the classical action Sci- Since the 
equation of motion is trivial, the classical solution satisfying our chosen boundary conditions 
is easily found, 

x(t) = x f + (t-t f f° + * { * f -*°\ (5.16) 



i + irzitf - u) 
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where the subscript / denotes final quantities. The classical action (15.91) . with (I5.10p . is 
then 

Sci(Xf ' tf) - Tk i + &(t f - ti ) ■ (5 - 17) 

Substituting this expression in (15.151) . we find 

A(t) = JL- (5.18) 

V 1 + 

with N a constant determined by normalising the wavefunction, N = e~ v ^ Ll '/ h2 / (2irL 2 )^ . 
Equation (15.61) then gives the familiar free-particle solution to the time-dependent Schrodinger 
equation, describing a spreading wavepacket. 

There are several points worth noting. First, the relevant classical solution is generally 
complex, and only real in the special case where Xf lies precisely on the classical trajectory, 
Xf = x c + p c (tf — U). Second, in general the solution does not start at x c , nor is its 
momentum p c . The initial value of x at t — t\ is 

x f - x c -p c (t f - U) 
Xi = x c + — 7—— — — — J -, (5.19) 

and value of the conserved momentum is 

If Xf is chosen to lie ahead of (or behind) the classical solution then the imaginary part of 
Xi is negative (or positive) and the complex classical trajectory lies in the lower (or upper) 
half of the complex x-plane. The initial condition ( 15.11) implies that pi — p c is directed at 
right angles to Xi — x c in the complex x-plane. Finally, the magnitude of the imaginary 
part of set by the difference between Xf and the value predicted by the classical 

solution, Xf — x c — p c (tf — U), times the dimensionless number (h/L 2 )(tf — ti) which is the 
fractional spreading of the wavepacket, AL/L ~ Ap{tf — U)/L where Ap ~ h/L is the 
initial spread in momentum implied by the uncertainty relation. If the fractional spreading 
of the wavepacket is small over the relevant interval of time, then for typical values of Xf 
the imaginary part of the relevant classical solution will also be small. 

We can easily generalize to a free particle confined to the positive half-line by a brick 
wall located at x = 0, forcing the wavefunction to vanish there. We can satisfy this 
boundary condition with free particle wavefunctions by the method of images. Assume that 
the particle starts out described, to a good approximation, by a Gaussian wavefunction 
obeying (15. ip . with x c large and positive and p c negative. To ensure the wavefunction 
vanishes on the wall, we subtract the mirror image wavefunction obtained by reflecting the 
initial conditions, x c — > —x c , p c — > — p c . Squaring the complete wavefunction, one finds the 
probability density 

'2xp c L 2 



ty\ (x } t) = C [e ^ +e ^ -2e ^ cos 



ha 1 



(5.21) 
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with C = 1/V27TCT 2 , cr 2 = L 2 + h the growing dispersion of the Gaussian and p c (t) = 

p c — x c h 2 (t — tj)/(4L 4 ). We minimize the spread of the wavefunction over the time t — t+ of 
interest if we choose L = \Jh(t — ti)/2. 

At early times, the probability density (15.211) is well-approximated by the first term 
alone, namely a slowly spreading Gaussian centered on the first classical trajectory. When 
the wavepacket hits the wall, at t — t j ~ —x c /p c , interference between the incident and 
reflected waves leads to strongly oscillatory behavior in both space and time. It is this 
regime which will be most closely analogous to the bounce from the singularity in our 
problem. As time runs forward, the second classical solution, obeying the mirror image 
initial conditions, now determines the dominant piece of the wavefunction. The probability 
once again assumes a smooth Gaussian form but now moving away from the wall. Similar 
behavior will be encountered below, when we deal with "reflection from infinity" in our 
problem. 

The semiclassical formalism we are using emphasizes the similarities between the quan- 
tum and classical problems, but their differences emerge in two interesting and important 
ways. First, more than one classical solution can be relevant to determining the wavefunc- 
tion. In the brick wall problem, two are needed and it is the interference between them 
which forces the wavefunction to vanish at the wall. Second, the relevant classical solutions 
are generically complex, and, in the brick wall problem, generically avoid the location of 
the wall in the complex x-plane. 



5.2 The Self-Adjoint Extension via Complex Classical Solutions 

In our problem, the homogeneous mode of the scalar field is described by the action for a 
quantum mechanical particle, 




(5.22) 



in which V(x) satisfies two general properties. First, it is unbounded below and falls faster 
than quadratically at large x so generic real trajectories lead to finite-time singularities. Sec- 

3 

ond, V (x)/ (— V{x))' 2 tends to zero at large x so the semiclassical approximation becomes 
exact there. As discussed in Subsection 14.11 under these two conditions the singularity 
can be quantum mechanically resolved by restricting the allowed stationary states to those 
proportional to a fixed real function of x at large x. For example, for V = — \\x A the 
allowed energy eigenfunctions behave as 

x — ► oo, (5.23) 

with a an arbitrary constant phase labeling the particular self-adjoint extension. 

For the case we are interested in, the use of the energy eigenfunction basis turns out to 
be an inessential complication. We are interested in the time evolution of a non-stationary 
state in which the homogeneous component of the field <fi is localized around some small 



ip%(x) ~ x~ 



cos 




2 3h 



+ oi 
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value in propagating the state across the quantum "bounce." In principle, we could compute 
the quantum evolution by expressing the initial state in terms of energy eigenf unctions, but 
it is far more convenient to describe the process using Gaussian wavepackets. As we shall 
now show, the self-adjoint extension may be implemented directly using complex classical 
solutions, for a class of potentials including those of interest here. 

The first step is to represent the theory (15.221) in the complex x— plane, viewed as a 
Riemann sphere including the point at infinity. The simplest potentials to deal with are 
those which are meromorphic in the region containing the singularity and symmetric under 
reflection through the singularity, x — > —x. The potential V(x) = — \\x A satisfies these 
requirements and we shall start by considering this case. 

Consider an initial Gaussian wavepacket with coordinates (x c ,p c )and spread Ax = L. 
The initial condition (15. ip provides one complex equation for x^ and = A second 
equation is obtained from energy conservation: since e = \x 2 — \\x Al is constant, we have 



For each value of Xf, and at fixed tf, equations (15. ip and (15.241) provide two complex 
equations for the unknowns Xi and Xi (or equivalently, e), which completely specify the first 
classical solution Xi(t). 

By analogy with the brick wall problem, we construct a second, "image," classical 
solution X2{t) by reflecting the initial condition (15. ip through infinity, 



but imposing the same final condition X2j = %ij — Xf at time tf. As for the brick wall, we 
expect the corresponding semiclassical wavefunction to start out "behind" the singularity, 
moving towards it, and to emerge from it. 

Now we want to show that the full semiclassical wavefunction for the self-adjoint exten- 
sion is just the sum of two terms 



each determined by the appropriate classical solution. As in the brick wall problem, the 
boundary conditions are satisfied by the method of images, due to the interference between 
the two terms. The main difference here is that whereas we had to take a = tt/2 to ensure 
the wavefunction vanished at the brick wall, in the present problem we can take any real 
< a < 7r and still ensure the Hamiltonian is self-adjoint. 

To show that the Hamiltonian is self-adjoint, we simply compute the wavefunction (I5.26P 
in the limit of large Xf. In this limit, the second classical solution becomes the reflected 
image of the first, X2(t) = —Xi(t), and this leads to the desired asymptotic behavior. 
Explicitly, we express the action in canonical form, 




(5.24) 




(5.25) 




(5.26) 




(5.27) 
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Because the two solutions are related by simple reflection in the limit of large Xf, they 
involve opposite momenta in this limit: in the first solution px = -K 2{e\ + tAx 4 ), whereas 



p 2 = — y2(e2 + j^x 4 ). Substituting into (15.271) and integrating by parts to extract the 
leading dependence for large Xf, we find 



Si - \l^ + j{tf-U) + 0{x?) 

3 




S 2 -> -^^L + ^tf-t^ + Oix] 1 ), xf^oc (5.28) 

where eo is the limiting energy of either solution as Xf tends to infinity, and = 
eo + O^xJ 1 ), as can be straightforwardly shown from equation (15.241) . The semiclassical 
wavefunction ( 15.261) thus behaves as 

* a (x, t) ~ e -^ e qe t/n x -i cog i ^ ^ z_ + a j ) x oo. (5.29) 

Since this is proportional to (15.231) . it is obvious that ^> a -^ip% — il)%-^^! a vanishes as x — > oo, 
for any energy eigenfunction Thus our time-dependent wavefunction \l/ a (x, t) lies within 
the domain of the self-adjoint extension. 

In this section, we have dealt only with the homogeneous mode. The analysis of Subsec- 
tion 14.21 shows how this can be extended to include inhomogeneous modes, in the presence 
of a UV cutoff. By changing coordinates to the (r, x s ) coordinates described there, one 
can show that the asymptotic behavior of the classical action near the singularity locally 
takes the form of (15.281) and hence that the self-adjoint extension is implemented point by 
point in the nonlinear theory, in the presence of a UV cutoff. Our calculation of the final 
wavefunction \l/(t/, 4>, <50(x)) will, for values of <fi carrying most of the probability, turn out 
to be insensitive to the UV cutoff. Hence we can describe the most probable final states 
in a manner independent of the UV cutoff, although we shall not be able to do so for all 
values of d). 



5.3 Dealing with Branch Cuts in V(x) at Complex x 

The above discussion is straightforwardly generalised to potentials whose analytic continu- 
ations are meromorphic in the region of the complex plane containing the singularity and 
which are symmetric under reflection through it. Some additional considerations are needed 
to deal with potentials, like that of interest in this paper, which possess branch points in 
the relevant region of the complex plane. 

Let us start by considering the general problem of constructing the semiclassical wave- 
function in a theory where V{x) possesses branch points somewhere off the real axis (Figure 
[6]). Since V(x) is real on the real x-axis, the branch points will occur in complex conjugate 
pairs. The branch points imply that V(x) is only single valued on a nontrivial cover of 
the complex x-plane consisting of a set of Riemann sheets glued together across the branch 
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Figure 6: Two complex classical trajectories passing on either side of a branch point in 
V(x), somewhere in the complex x-plane. 

cuts, which indicate the passage from one Riemann sheet to the next. A classical trajectory 
which runs into a branch cut should be continued onto the next Riemann sheet, upon which 
the potential V(x) takes a different set of values at the corresponding complex coordinate 
x. For example, if V(x) = ln(x 2 + a 2 ) with a real, then there are two branch points at 
x = ±ia. If we cross a branch cut emanating from one of them the potential changes by 
±27u as we pass onto the next Riemann sheet. 

We are interested in potentials V(x) which are real on the real x-axis where we want 
to compute the wavefunction. However, if V(x) has branch cuts then it is possible for a 
classical trajectory to cross a branch cut and go to the next Riemann sheet, on which the 
potential will no longer be real on the real x-axis. In some cases the branch cut can moved 
out of the way and then we can simply ignore it. However, in general this will not be 
possible. For a given initial condition, as we vary the final coordinate Xf or the final time 
tf, at some critical value the relevant classical solution may encounter the branch point. 
Computing the wavefunction at values of Xf or tf just beyond the critical value then involves 
dealing with solutions which cross the branch cut. If we continue the solutions onto the 
next Riemann sheet so that they again end on the real x-axis, the resulting semiclassical 
wavefunction will solve the Schrodinger equation with the wrong potential V(x) - the one 
appropriate to the next Riemann sheet. Instead, we have to change the initial condition of 
the relevant classical solution so that it begins on the next Riemann sheet and ends on the 
correct Riemann sheet, where V(x) takes the physical value on the real x axis. 

Luckily, there is a simple procedure which allows us to calculate the correct initial 
condition for classical solutions which cross a branch cut. It relies on the fact that as long 
as V(x) is sufficiently well-behaved and, in particular, has no 5-function contributions, the 
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Schrddinger wavefunction \1/ and its first space derivative (d/dx)^f(t,x) must be continuous 
functions of x. Recalling that within the semiclassical approximation, (h/i)(d/dx)^(t, x) ~ 
p(x)^(t,x) at leading order, we see that the momentum of the classical solutions being 
considered must be continuous as the final coordinate x is varied. 

For a given initial condition (15. ip . with fixed x c and p c , consider decreasing Xf along 
the real x-axis so that at some critical value x*^ the classical solution hits the branch point. 
First consider a classical solution ending just before the critical value, at x*j — e say, with 
e infinitesimal and positive. This solution passes just below the branch point, and ends 
on the real x axis with some final (generically complex) canonical momentum p*j. Now 
consider a solution ending just beyond the critical point, at x*j + e. If the semiclassical 
wavefunction ty(t,x) is to be continuous and have a continuous first derivative at x**, then 
in the limit of small e the classical solution ending at x*j + e must have the same final 
momentum p*j. With its final coordinate and momentum now determined, we can run 
this supercritical solution back in time from tf to ij. It will pass just above the branch 
point, across the branch cut, and onto the next Riemann sheet. The potential V(x) will be 
different on this sheet and hence the trajectory followed will differ from that of the original 
solution. Likewise, the classical action calculated for this second solution will differ from 
that of the first. Running the second solution back to U, we find the initial values x, and 
Pi and hence the initial value of the coherent state parameter in ( 15.11) . This new value is 
to be used as the initial condition for the classical solutions determining the wavefunction 
beyond x**, assuming these solutions cross no additional branch cuts. Besides shifting the 
initial condition for these solutions by a discrete amount, we must also add a (complex) 
constant to the action S in order to ensure that the wavefunction is continuous across x*j. 
This can be thought of as a "boundary term" associated with passing around the branch 
point and onto the next Riemann sheet. Finally, although we have described crossing the 
branch point in terms of varying the final position Xf, we could have just as easily varied 
the final time tf. The same procedure we have given then ensures that the wavefunction is 
continuous when considered as a function of time. 

To summarize, when the classical solutions of interest move across a branch point, we 
have to make a one-off change to both their initial condition (coherent state parameter) 
and their action. Both changes are determined by the continuity of the wavefunction and 
its first x derivative on the real x-axis. 

5.4 Resolving the Branch Points arising from Renormalization 

For the theory at hand, we have a double problem. The renormalized effective potential 
for the scalar field <p, when written in terms of the canonically normalized variable x, is 
proportional to — x A / ln(Ca;). We can examine the behavior at infinite x by setting \ = x~ l , 
so the potential takes the form x~ 4 / ^ n {x/C), with C a constant. The logarithm gives rise 
to a branch point at x = which vitiates the construction of the self-adjoint extension 
given in Subsection 15.21 

In order to apply the method of images, the analytically continued potential must be 
even under x ~^ ~X- This will be the case if V(x) is defined with two branch cuts, running 
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along the positive and negative portions of the real %-axis respectively. The logarithm may 
then be taken to be real just above the positive real axis and just below the negative real 
axis. However, to solve the correct Schrodinger equation we must ensure that the potential 
V(x) is real at the endpoint of both classical solutions - one of which hits the positive real 
X-axis from above, and one which hits it from below. Unfortunately, we cannot have it 
both ways - we cannot simultaneously symmetrize the potential under x ~~ > ~X and have 
it real on both sides of the positive real axis. 

The solution of the problem is to separate the branch cuts from the fourth order pole. In 
Appendices O and Owe show this is naturally accomplished by renormalizing the theory at 
a large but finite cutoff p, constructing the self-adjoint extension from the relevant classical 
solutions and then taking the cutoff to infinity at the end of the calculation. Cutting off 
the high momentum modes turns off the running of the coupling A^ so that, for |0| ^> p, A^ 
tends to a constant. Therefore V(<j)) only has a fourth order pole and no branch point at 
<p = oo. The branch points associated with the logarithmic running of the coupling with \(f>\ 
have therefore been moved away from x — 0. In Appendix Owe study the running coupling 
in detail, showing that for finite p there are four branch points arranged in a symmetrical 
pattern roughly \Aoo/ P from x — 0, where Aoo is the value of the running coupling at 
= oo. Furthermore, the running coupling is finite and roughly equal to A m at the four 
branch points. 

We can now discuss the construction of the semiclassical wavefunction using two classical 
trajectories as in Subsection 15.21 The branch cuts in V can be placed in such a way as to 
respect the symmetry x ~ * ~X while leaving V real on the real x axis. As we shall see 
in Subsection 15.61 for very short times and small Xf both classical solutions run close to 
the real axis in the small x region . Both miss the branch cuts and the method of images 
works just as in Subsection 15.21 Now consider increasing the time tf — U or the value of Xf, 
so that one of the two solutions encounters a branch point. As explained in the previous 
subsection, we must correct the initial conditions and the action for all subsequent solutions 
in order to ensure continuity of the wavefunction and its first x derivative. This procedure 
may be complicated to implement in practice, but its existence guarantees that, at finite 
cutoff p, we can construct a semiclassical wavefunction satisfying the correct Schrodinger 
equation and lying within the domain of the self-adjoint extension. 

We now want to consider the limit as p tends to infinity. We would like to know whether 
the complex classical solutions we need are sensitive to the details of the UV cutoff p . The 
fact that the effective potential is finite at the branch points and varies smoothly away 
from them on a given Riemann sheet suggests that the infinite p limit may be smoothly 
taken. Consider two neighbouring classical solutions, as in Figure El running back from 
neighbouring values of \f and possessing the same energy and final momentum. They pass 
respectively just above and below the branch point, with one of them crossing the branch 
cut. Their trajectories are each given by x = ^2 (e — V(x)), so the velocity difference 
grows smoothly away from the branch point, and the two trajectories gradually diverge. 
Since the difference between the trajectories builds up slowly, it is plausible that detailed 
structure of the potential near the branch points does not matter, and that the p — > 00 
limit can be taken before calculating the two relevant complex solutions. 
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In the limit of infinite p, the effective potential is proportional to x~ 4 /hi(x/C)- In 
order to respect the \ — > — x symmetry we place two branch cuts along the positive and 
negative real %-axis respectively. Defining ln(z) to be real for positive real z and to possess 
a branch cut along the negative real z-axis, h(x) — \ ( m (x) + ln(— x) + in) is even and 
real just above the positive real x _ax is or just below the negative real x _ax is. Likewise, 
h(x) = h(x) — in, which is the extension of l\ to the next Riemann sheet, is real just below 
the positive real axis and just above the negative real axis. 

Consider the physical situation we are interested in. We start the field in a Gaussian 
wavepacket rolling gently downhill, at nearly zero potential energy. We would like to 
calculate the wavefunction after one bounce, when the field has returned to a state of nearly 
zero potential energy. The self-adjoint wavefunction comprises two pieces, corresponding 
to the first and second, "image," classical solution respectivly. At the times of interest to 
us, the latter part of the wavefunction dominates, depending on the classical solution which 
started "behind" the singularity and passed around x = before ending at large, positive 
Xf- (Figure M below shows such a trajectory.) In this case, we should define the effective 
potential V using h(x) to ensure it is real as x approaches Xf from below the positive 
real axis. This second classical solution is defined with a certain coherent state parameter 
defining x c and p c appearing in the initial condition (15.11) . However, this is not just minus 
the coherent state parameter of the first classical solution, because the two solutions evolve 
on different Riemann sheets, where the logarithm takes the values defined by l\ and 1% 
respectively. The coherent state parameter of the first classical solution is the one which is 
most directly specified by the initial conditions. 

In order to relate the two coherent state parameters, we must study the transition from a 
regime in which both solutions evolve on a single Riemann sheet. This is pictured in Figure 
[7J Consider the wavefunction at a time £/ — 1< given by roughly twice the classical roll-down 
time and focus on the behavior as we vary the final value Xf- As X/ is lowered, the second 
classical solution deforms smoothly to the lower curve shown in Figure UJ This solution has 
to reach small X/ in roughly twice the normal classical time. It does so by "loitering" in the 
lower half x-pl an e with its real part first crossing zero and "rolling uphill" before turning 
around and rolling back to small x- Recalling the location of the branch points in the finite 
p theory, we can draw the branch cuts so that the potential is symmetrical and the first 
solution is just the reflection of the first, Xi(t) — ~X2(t), as long as Xf is sufficiently small. 
It is clear that this entire family of classical trajectories never encounters a branch point. 
So we can use the potential with the log defined by h(x) f° r an °f them. 

Now let us study the first classical solution. In the limit of small Xf an d in the finite p 
theory, we can draw the branch cuts as shown in Figure[7Jso that the potential is symmetrical 
and the first solution it is just the reflection of the first, Xi(£) — ~~ X2(t)- As Xf is increased, 
the first solution crosses the branch point and runs onto the next Riemann sheet. Increasing 
Xf further, the first solution has crossed both branch points. Now consider taking the 
p — > oo limit but still focusing on the values of Xf ne ar zero, for which the solution passes 
just above or below the branch points. The solution passing just below the branch points 
should, in this limit, evolve for almost the whole of its length on the sheet where the 
logarithm is Whereas the solution passing just above both branch points should evolve 
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Figure 7: The two classical solutions used to determine the semiclassical wavefunction are 
shown, where the effective potential is calculated with a large but finite cutoff p so that it 
is regular at x = but has four branch points arranged symmetrically around the origin. 
First, the two solutions are shown at final times tf corresponding to the classical time for 
a full bounce, but for Xf being close to zero. These solutions have a large imaginary part, 
and circumnavigate the branch cuts. As is increased, the first solution encounters the 
branch cut, and \i m ust be modified as explained in the main text. 
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for almost the whole of its length on the sheet where the logarithm is given by l±. Since 
the logarithm is very large for small |x|, the effective potential on the two sheets hardly 
differs at small \x\- The difference between the two classical trajectories builds up along 
the solution as they run to larger values of |%|. 

In Appendices O and [TJwe justify the neglect of all of the details of the finite-p potential 
in the vicinity of the branch points by first showing that all of the features of the potential, 
including the branch cuts, are finite and suppressed by a positive power of the perturbative 
coupling constant. Then we show that such small features of the potential, only relevant in 
the small x region of the trajectory, have a negligible effect on the trajectory. This justifies 
our direct use of the infinite-p potential in order to compute the classical solutions. 

We calculate the coherent state parameter for the first solution as follows, using the 
infinite p limit directly. We first find the solution which, with initial conditions given by 
reflecting those of the second classical solution and with the I2 definition of the logartithm, 
hits X/ = at the time t/ of interest. We then measure \ and \ at some time very close 
to tf on this solution. Then, taking these same values we switch to the Zi definition of the 
logarithm and run the solution backwards in time to £j. The values of \ and \ measured 
there define the coherent state parameter for the first classical solution, and this is used to 
define the first piece of the semiclassical wavefunction for all positive Xf- 

In subsequent subsections, we will implement this procedure in some detail. First, we 
will find the relevant complex solutions analytically, in the approximation where the loga- 
rithmic dependence of the potential can be ignored. Then, after explaining how the nonzero 
momentum modes can be simultaneously described using the semiclassical approximation, 
we shall move on to a full numerical calculation of the relevant complex solutions. 

5.5 Complex Solutions at Fixed Coupling 

In this subsection, we want to explicitly construct the complex classical solution describing 
a bounce off the singularity at x = 00. We shall in the first instance ignore the logarithmic 
running of the quartic coupling A^,, treating it as a constant. After a field and coupling 
redefinition to remove factors of the spatial volume, the action is just 




(5.30) 



Furthermore, we can remove A by setting x — a/2/A 



x 



-1 



so 




(5.31) 



We want to find classical solutions satisfying (15. ip . which reads 



o ■ r 2 , 1 r,- t2Xc . 1 , 

- 2iL —2 H = —2tL — H , t 



Xi Xi X c Xc 



t = u 



(5.32) 




(5.33) 
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with e a constant. Integrating, we obtain 




y/l + 2e X 4 ' 



dx 



(5.34) 



Note that L is the width of the Gaussian in the variable x: the width in the homogeneous 
component of the field, 0, is A0 = Lj\/V^. Eliminating Xt i n favour of e using (15.331) . 
equations (15.341) and (15.321) provide two complex equations for two complex unknowns, x% 
and e. 

We shall solve the equations in the regime most relevant to the description of the final 
state after a bounce. For this, we need the second, "image" solution, with initial conditions 
that Xc is large and negative and Xc ~ 1 so that we are close to the zero-energy scaling 
solution. We want to evolve to a final x / which is large and positive so that we have crossed 
the singularity and returned to the configuration where the homogeneous component is 
close to its initial position. 

We shall look for complex classical solutions close to the real classical solution x — t, 
hence the imaginary part of x will be small. Setting x% = r % + ^ s ii we first calculate the real 
and imaginary parts of (15.321) . To linear order in s, these read 



where r c = Xc (which is real). Since the energy of the solutions we are interested in is small, 
it is reasonable to expand the square root in (I5.34j) in ex 4 , assumed small. To leading order 
in e we obtain 



Now we make several further approximations, again related to our assumption that the 
complex solutions of interest are close to the real classical solution. First, we assume s is 
small throughout, so we need only work to first order in s. Second, we assume that x% is 
close to Xc so that = x c (l — 5), with (5< 1. And finally, we assume that since 8, e and 
s are all small, we can ignore any term involving products of them. 
The initial conditions (I5.35P now give 




(5.35) 



tf-U ~x/-Xi- ^{x)-x\)- 



(5.36) 




(5.37) 



where e = e + iei, with e and e\ real. In the same approximation, (15.361) reads 




(5.38) 
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Equations (15.371) and (15.381) provide four equations for four unknowns, s, 6, eo and e\. 
We consider a final Xf close to — Xc-, so we are looking at the final wavefunction around the 
initial value of <ft, to first order in small quantities. We find 

5L 2 J 1 



h Xcl + L\h Xc )- r 
L 2 



5 
5 

ei = e rSj, (5.39) 

X f ~ Xl 

where J — tf — t% — (Xf ~ Xc) measures the deviation of the argument of the wavefunction 
from the classical zero-energy solution. 

Of particular interest to us later will be the minimum value of |%| attained, since that will 
act as a UV cutoff for the quantum production of inhomogeneous fluctuations. Evaluating 
(15.341) . with ti replaced by t, we find s is quintic in r and 

Smin = -~^X) = 7}Si, (5.40) 

with Si given by f!5.39[) . Neglecting the quantum spreading term, we have 

Smin ~ — — L 2 , (5-41) 

2 Xc 

where, as we recall, L is the width of the initial Gaussian wavepacket in x, Xc is the initial 
value of the homogeneous mode of x an d J is the deviation of Xf horn the final value 
expected from classical evolution. 

Let us now translate this result into the situation of interest. We want to consider an 
initial Gaussian wavepacket with nearly zero energy, but over the barrier, i.e., centered 
on 0c ~ R-Ads- We parameterize the initial width A0 of the Gaussian in terms of 
the width A0 m j„ which minimizes the spread of the wavepacket over the entire duration 
of the putative bounce, a time of order RacLS- It makes sense to perform the calculation 
in terms of the variable Xi defined so that x = \A is the classical zero-energy scaling 
solution. In the presence of the running logarithm we thus define x = d<p{2/ X^) 1 / 2 ^ 2 . 
The final variance in x is then approximately given by Axf + R 2 AdS {Axi) 2 ■ Here Axi ~ 
(A0)A; V V 2 = (Att^^A^V 2 ~ (A^-^A^V 2 ~ ^R/(A Xi ). Thus the final 
variance is Ax 2 + X 2 R\ dS /{Axi) 2 , minimized by Axi ~ Rms < Xi ~ RacIS- The 
corresponding value of A(p min is R^d S <C 4> c - 

We shall write the width of our initial wavepacket A0 = WA<f) m i n , where W is a 
dimensionless parameter. From (15.411) . the value of s m i n is, for typical values of Xf where 
the wavefunction has most of its support, of order W 2 RAdS^Xi/Xi ~ W 3 Xl/ 2 RAds- As 
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we shall see later, this small imaginary part for the classical homogeneous background 4> 
provides a UV cutoff on particle production across the bounce, for a range of final values 
of <pf carrying most of the probability. 

Before we calculate particle production in this background, however, we want to study 
the large-0 behavior of the wavefunction at very short times, in order to justify our claim 
that the quantum spreading cannot be ignored, in principle, for arbitrarily short times 
and hence the dynamical evolution cannot even be discussed without imposing a unitary 
boundary condition at = oo. 

5.6 Behavior of the Wavefunction at Large <p and Small t 

One of the unusual features of the system we are considering is that even though the 

quantum evolution of any state is perfectly well-defined, neither the expectation value of 

the homogeneous mode <p, nor any of its correlators, exist. Hence the usual description of 

quantum field theory in terms of correlation functions of fields, does not, strictly speaking, 

make sense. Instead, one has to deal directly with the wavefunction itself. That the 

correlation functions are in general ill-defined may be seen from (14.61) . for example. Even 

though the energy eigenf unctions are all normalizable since |\l/| 2 ~ x~ 2 cos 2 (Ax 3 + a) at 

i 

large x, the expectation of x — <pV 3 2 is a meaningless logarithmically divergent integral 
~ j (dx/x) cos 2 (Ax 3 + a) in any one of them. 

At first sight, this behavior appears puzzling. There is no problem in setting up a 
Gaussian wavepacket at time t = 0. In such a wavepacket, the expectation value of x or 
any integral power of x is well-defined. Likewise, the time derivative (d/dt)(x) = (x) = (p), 
is also perfectly well-defined. Using the Heisenberg equations of motion, one can express 
the n'th time derivative of (x) as an expectation value of a polynomial in x and p, whose 
value is perfectly well-defined at the initial time. Hence, formally, (x) is a function whose 
value and time derivatives to any order are finite at t = 0. Nevertheless, (x) does not 
actually exist, for any t > 0. 

We can understand what happens by calculating the semiclassical wavefunction using 
complex classical solutions. Consider the case where we are given a Gaussian wavepacket 
for x whose initial coordinate x c is small, and initial speed x c is given by the zero-energy 
scaling solution. We want to study the form of the wavefunction at very large Xf, at small 
positive final times tf. (We set ti = in this subsection.) 

As in the previous subsection, we change variables to x — (2/ A) 1 / 2 ^™ 1 , for which the 
initial conditions are given by (15.321) with ^ c C 1 and Xc = — 1- Again, conservation of 

energy yields 

X = - V / l + ex 4 , (5-42) 

with e a constant. The sign is chosen because the solution we are interested in runs from 
some large initial x to a very small real final Xfi after a short time tf. The form of the 
solution is illustrated in a specific numerical example in Figure [H 

Recall that the classical solution does not need to start at Xc- If we take Xc large, 
the second term on the right hand side of (15.321) dominates. We now solve for the initial 
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Figure 8: Complex classical solution describing the large-0, small time behaviour of the 
wavefunction for the homogeneous background. The trajectory of the variable \ is plotted 
in the complex x-plane. 
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(5.43) 



which is small and nearly negative imaginary if Xi is nearly real. In order to achieve 
consistency with (I5.42p . we must take Xi close to the branch point in the right hand side 
of f)5.42p . Xi ~ Xb = (— e) _ 4. If e = eo + ie± with eo large and negative, and e± small and 
positive, then xb ~ (— eo) _3 (l — ^i/(4e )) has a small positive imaginary part. Setting 
X = Xb + Sx, (|5-42p reads 



X = -2(-5 X )H-eY 



(5.44) 



Thus, if x starts out just to the right of xb, *-e., with 5x small and positive, then from 
(15.441) . Xi is negative imaginary as required. As x runs downwards towards the real x-axis, 
from (15.441) . x turns towards the origin so that x heads to small, real values. Solving (I5.43P 
and (15.441) for Sx, we find 

x * = - Tek) ■ < 5 - 45 > 

just to the right of the branch point. 

If the real part of the energy eo is large and negative, as we have assumed, then Xi is 
small. We are assuming that Xf is even smaller. But then ( 15.341) becomes 



tf - U 



dx 



o V 1 - WxJ 



Xi 



dz 



o VT^ 



kx 



(5.46) 



where k = \A"(r(|)/r(|) = 1.31103. . . is the complete elliptic integral of the first kind 
with parameter %. From this value of Xh (15.431) gives Xi- Neglecting terms down by Xi/Xc 
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and using (I5.46P we obtain 



*=-^!=(^_ (5 . 47) 

at the initial time. Armed with this classical solution, we can compute the classical action, 
with the boundary term (I5.10p . The integral is dominated by times t near t*, the time 
when the solution, if continued past tf, would encounter the singularity: 

r** 1 iL 2 x 2 

Sc, = l dt W^W + ^r- (5 ' 48) 

The integral is re-expressed in terms of Xf using t* — tf + (2/X) l l 2 x^ 1 ; the boundary term 
is computed from (15.471) . We obtain 



1/2^3 ^2 

J7 



iSci (A/2) IZZ _^ 

h 3h 2h 2 \L 2 t 2 f ' 1 > 



Now we can read off the behavior of the wavefunction at very large Xf for small final 
times tf. Including the contribution from the second classical solution, where the initial 
conditions were reflected through x = 0, we find for the total wavefunction at large x and 
small times tf, 

*(x, t f ) ~ x^exp [~ 2h2 K XL 2 t 2 ) cos \\l 77^ + a ) • ( 5 - 5 °) 

_ .—2 

Note that the exponential factor ~ e / , with c a constant, is indeed a function whose 
value and time derivatives of arbitrary order vanish at t = 0. Nevertheless, the function is 
nonzero for any nonzero t. Thus, although correlators of x or p are well defined at t = 0, 
for any positive t the probability acquires a tail at large x proportional to x~ 2 times an 
oscillatory factor which averages to |. Thus for any positive time t, the expectation value 
of x n with n > 1 does not exist. 




5.7 Complex Solutions for <p with a Running Coupling 

As we shall see later, a critical factor in our calculations will be the running of the quartic 
coupling A: although this is a weak effect, without it, the particle production would vanish 
at linearized order. We start from the effective potential obtained in (13.221) . 

S = j d*x (~(<90) 2 + ^0 4 ), (5.51) 

and we again ignore the curvature of the S* 3 because this is unimportant near the singularity. 
We are interested in the regime where the logarithm I in the denominator of A^ is large: in 
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this regime, the logarithm varies slowly with <p and we can solve the equations of motion 
in an expansion in inverse powers of I. 

As before, we first change variables to x — (2/Ao) obtaining the action 

5=f [^(-(dxf + r 1 ). (5.52) 
Ao J X 

In order to minimize the effect of the logarithm, we redefine the field variable, setting 

d\i = I* (x) d>X, ( 5 -53) 
where I = In ((2/X )^/(NMx)). This is solved by 

Xi = X ; HI) = 1 + I- -L + A.... (5.54) 

We obtain the following action for xii 

If d A xlh 4 (l) , x9 v . 
S=t- {-(dXi) 2 + 1) , (5-55) 



Ao J Xi 

in which form it is clear that at small x, one obtains a large classical action and therefore one 
can expect the semiclassical approximation to become accurate. The equation of motion 
for xi is 

- X id 2 Xi + 2{d X if + 2 = -2g(l) ((d X i) 2 + l) , (5.56) 

where g(l) = h{l) - 1 - \l~ l h{l) = - f/" 1 + \l~ 2 + ...). 

From (I5.56P one observes that xi = ±(^s — t) with constant t s still solves the equation, 
just as before. By redefining the field as we have, we have removed the logarithmic cor- 
rections from the background solution. We could proceed to write the general solution in 
the vicinity of a spacelike singularity by simply modifying (14.171) so that each coefficient of 
T n , n > 2 and r n lnr, n > 5, becomes a series in inverse powers of I. Near the singularity, / 
diverges so these logarithmic corrections become less and less significant. 

Figure [9] shows an example of a complex solution to equation (15.561) describing the 
homogeneous mode Xl(t)- This solution corresponds to the "image" solution which starts 
out behind the singularity at x = and runs to a positive real final value Xf- Notice 
that the behavior of the complex solution is very simple in the vicinity of the origin: xi is 
well- approximated in this region as xi = t — it with e a constant. In the next section, we 
shall discuss the evolution of the inhomogeneous modes in this background. 



6 Quantum Evolution of the Inhomogeneous Modes 

Having discussed the quantum description of the homogeneous mode of the scalar field, we 
turn to the quantum evolution of the inhomogeneous modes, 

0(t,x) = ^0 k (t)F k (x), (6.1) 

k 
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Figure 9: Complex classical background solution in which particle production is computed, 
including the effect of the running coupling. The parameters used were NM = 10~ 5 , L = 1, 
U = -10, n = 10" 5 , Si = -0.4, tf = 8.4. 

where lk(x) are the spherical harmonics on S 3 , and k are the corresponding eigenvalues. 
We shall be interested in modes whose wavelength is much shorter than the radius of the 
S 3 , i.e., RAds, so that we can approximate the spherical harmonics with plane waves e ikx . 
As explained in Appendix [Bj our calculational procedure is to integrate out the oscillatory 
modes, i.e., those for which \kt\ > 1, and to study the long-wavelength modes, i.e., those 
for which \kt\ < 1, in the linearized approximation where we neglect their backreaction on 
the homogeneous mode. As we shall see, particle creation will turn out to be finite and 
suppressed by powers of the coupling A^ (see ( I3.23P ). so this approximation is reasonable 
for small A^. 

We want to describe the homogeneous mode, and the k ^ modes, in a single semi- 
classical wavefunction, obtained from complex, spatially inhomogeneous solutions to the 
classical equation of motion for the field. In the case where the quantum spreading of the 
wavefunction is negligible for the homogeneous background, this treatment will reproduce 
the usual description of quantum field theory in time-dependent backgrounds, where one 
uses the Heisenberg picture to compute Bogolubov coefficients, and then studies particle 
creation etc. We shall describe the same phenomena in the more general context of a quan- 
tum mechanical background, using the Schrodinger wavefunction, calculated from complex 
classical solutions. Before proceeding with our specific case, let us then briefly review how a 
quantum field in a time- dependent background can be described in the Schrodinger picture 
(see, for example, Ref. [57]). 

6.1 Particle Production Using Complex Classical Solutions 

The description of linearised perturbations around a time-dependent, spatially homogeneous 
background can always be framed (after a possible field redefinition to put the kinetic term 
into canonical form) in terms of a quadratic action 




(6.2) 
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where q is the coordinate of the particular linearized mode being considered (for example, 
the real or imaginary part of the 0k mode coefficients in (16.11) above) and oo 2 (t) is a time- 
dependent frequency squared. 

We would like to describe the situation where the inhomogeneous modes start out in 
the incoming adiabatic vacuum. That is, we assume that for the modes of interest, u 2 is 
positive, and uj/uo 2 is far smaller than unity, so that the mode behaves like a harmonic 
oscillator with a slowly-varying frequency. We also assume that every mode starts out in 
its adiabatic ground state. This is a special case of the formalism described above, for 
describing the time evolution of Gaussian wavepackets. In the present case, the initial 
width of the Gaussian is determined in terms of the frequency to. It is straightforward to 
determine the classical action, including the boundary term needed to describe the incoming 
adiabatic vacuum: in canonical Hamiltonian form, 



Sci = J dt [pq - ~(p 2 + uj 2 (t)q 2 ) 1 + (6.3) 

where Pi and u>(ti) are the initial momentum and frequency. In order to calculate the 
Schrodinger wavefunction for the mode, the correct initial and final conditions are: 

p = q = iu(ti)q, t = U\ <? = <?/, t = t f , (6.4) 

with qj real. The equation of motion is just q = —u> 2 q, and we assume that at U the 
evolution is adiabatic. Hence the two independent solutions are proportional to e +% ^ udt and 
e -ifudt reS p ec tively. Clearly, the initial condition (16.41) selects the former - the incoming 
negative frequency mode, which we shall denote Rl~\ while the final condition specifies its 
amplitude. The desired solution is then just 



We then easily compute the classical action (16. 3p . using integration by parts and the equa- 
tion of motion to obtain ^ ^ 

*~e^, S C i = l^fjq}. (6.6) 

Rin 

This formula agrees with that obtained by solving the time-dependent Schrodinger equation, 
given for example in Ref. [57]. 

Having obtained the complete qj dependence of the Schrodinger wavefunction, we would 
like to calculate the number of particles in the final state. The calculation is straightforward. 
First, we construct the creation and annihilation operators in the Schrodinger picture, at 
the final time tf\ 

ip f + u f q f 1 f d \ t -ip f + u f q f I ( d 

a = . — — . h- h LOtqt , a' = . — — — . — n— \- LOtqt 

^2uj f h y/2u f H V dq f 1 ) y/2u]h ^2uJ f h \ dq f 

(6.7) 
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Since ujj = u)(tf) is real, a' is the Hermitian conjugate of a, and the number of particles in 
this mode in the final state is 



This is our final result, giving the expected number of particles in mode k in terms of the 
behavior of the incoming negative frequency mode at late times. In the next subsections, 
we shall use this formula to derive the spectrum of produced particles. 

In passing, let us note the relation between (16.81) and the usual Bogolubov coefficients, 
calculated by evolving the inhomogeneous field modes in a real classical background. If 
the mode evolution is adiabatic at the final time tf, we can express the incoming negative 
frequency mode in terms of the outgoing negative and positive frequency modes, 



+if Ujdt 

p(-) . 1,1 

\/2uj 



ae 



+i /' Ludt _|_ jJ(D—i /' wdt 



2ui 

with uj nearly constant and real. Substituting (16. 9p into (16.81) we obtain 



*-t,, (6.9) 



l«l 2 - w 



If uj is always real, then the conservation of the Wronskian R^*R^~' — R^ 1 * R^ 1 yields 
|a| 2 — |/?| 2 = 1. However, when u is complex, we have to use the more general result (16.101) . 

In this subsection we have shown how complex solutions of the classical field equations 
may be used to compute the Schrodinger wavefunction for the fluctuation modes. In the 
following subsections we shall solve for the positive frequency mode and explicitly determine 
(n) for various particle species in the dual theory. 

Before doing so, let us make the following general point. When the Higgs field $i, whose 
magnitude we have parameterised by 0, rolls to large values, certain gauge, Higgs and Fermi 
particles acquire time- dependent masses M. of order g<p, where g is the Yang-Mills gauge 
coupling. A dimensionless measure of the slowness of increase of A4, i.e. the adiabaticity 
of the change in A4, as the singularity approaches is M./M 2 . If this quantity is small, one 
expects particle production of particles of mass M. to be exponentially suppressed. In our 
case, using the scaling solution ~ (\ ( p/2) 1 / 2 (j) 2 we have M./M 2 ~ \f\>/g, which tends 
to zero as the singularity approaches because is asymptotically free (see (15.511) ). The 
situation is further improved when one takes into account the dependence on N, arising 
from the gauge group SU (N) of the dual theory. To study the theory at large N, one keeps 
the 't Hooft coupling g t = g 2 N fixed, while A = 2a 2 /N 2 ln(0/iVM) (see (I53TD and (153311 ). 
Hence the parameter measuring the departure from adiabaticity in the massive fields is 
M/M 2 ~ l/(g t N\n((j)/NM)y/ 2 . 
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Therefore, when the logarithm and N are large, it should be a good approximation to 
integrate out the particles with masses of order gcj) and consider only the remaining degrees 
of freedom with tachyonic masses (the fluctuations 5(f) in the rolling field itself), with light 
masses, like Higgs modes acquiring their mass only from the deformation, or zero mass, like 
the massless gauge bosons. 

In the following subsection, we shall derive the equations of motion for the relevant 
quantum field in each case and show that, in an expansion in inverse powers of the logarithm 
/ = \n(<p/NM), to lowest order there is no particle creation in the vicinity of the singularity. 
Then, in the subsequent three subsections, we shall calculate the particle production in each 
species at the first nontrivial order in an expansion in inverse powers of I. 



6.2 Equations of Motion for Inhomogeneous Fluctuations 

In the following subsections, we shall compute the evolution of the inhomogeneous modes 
and the resulting particle production. We focus on those wavenumbers which depart from 
adiabatic oscillatory behavior (or, in more colloquial language, which "freeze out"), in the 
regime of large where the classical background is well- described by the zero-energy scaling 
solution. 

Let us start by recalling the zero energy real scaling solution, 

A A 2r7 2 

f<p\ \ = l = HMNM), A = ^. (6.11) 

Note that for the final times we are considering, roughly the time taken for the field to 
roll classically to infinity and bounce back, the quantum wavefunction is dominated by 
the second, "image" complex classical solution in which <fi emerges from "behind" the 
singularity. We can integrate (16.111) to obtain 





, i 




where the singularity has been placed at t — and the series in Z -1 is obtained by repeated 
integration by parts. From this we obtain 

It is also useful to calculate the time derivative of I, 

i=~(l + ^ + -..)■ (6-14) 

As we have explained at the end of Subsection 15.71 the complex solution which describes 
the homogeneous mode in the vicinity of the singularity is obtained by replacing t with 
t — ie in the above expressions. 
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Figure 10: The frequency squared u 2 (see equation (16. 15ft ) from the full complex numerical 
solution compared with the approximation xi — t — ie. The two sets of curves show the 
real part of u 2 and its approximation (solid curves, red and blue), and the imaginary part 
of uj 2 and its approximation (dashed curves, green and yellow). The curves corresponding 
to numerical and approximate solutions are indistinguishable in the figure. Parameters are 
as in Figure [91 

The equation of motion for the inhomogeneous perturbation modes is obtained by com- 
puting and substituting (16.131) . One finds 



where in obtaining the second line we have used the definition ( 15. 54ft of xu which is well- 
approximated by t — ie in the complex classical background of interest (see Subsection 15.71) . 
A check on the accuracy of this substitution is given in Figure EH where we plot the 
numerically calculated complex frequency squared, u 2 in (16.151) . versus the approximation 
obtained by setting xi — t—ie, with e a constant. Clearly, this is an excellent approximation 
over a wide range of times, and hence for a wide range of wavenumbers k. 

We wish to solve equation (16. 15ft as a series in At lowest order, ( 16. 15ft becomes 



5(j) 




8<j) 



= — uj 2 S4>, 



(6.15) 



5(f) 



6 



-5<f> - k 2 5<P, 



(6.16) 



(t - ie) 
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which is regular for real t. The solutions are (t — ie)z times Bessel functions of order 
v — 5/2. The two linearly independent solutions may be taken to be: 

,m , ~ ( 1 1\ sin/ct , ~/l 1 \ cos kt 

/«=cosA;t + — ^, / (2) = sin fct - - + ^. 6.17 

J V(H) 2 3J kt J \3 (kt) 2 J kt V ; 

where t = t — ie. The incoming negative frequency mode, z.e., the solution tending to e lkt 
at large negative times, is just — 3(f^ — if^)e~ kt . Evolving this to large positive times, we 
find there is no positive frequency component and, from (16.81) . the net particle production 
is zero. To discuss the next order in Z _1 we shall have to solve (16.151) to this order and we 
shall do so in the following subsection. 

Let us now work out the equation of motion for the light Higgs particles, i.e., those 
acquiring masses solely from the deformation. The latter contributes a time-dependent 
positive mass squared |A^Tr($ 2 ) = i^A^ 2 to each of the Higgs fields i — 2, . . . 6. In the 
scaling solution for 0, the linearized equations of motion for the light Higgs fields are 

^ = ~ 5(t ^ e)2 ( 1 + l ~ l + • • 0$; - k 2 ^, % = 2, ... 6. (6.18) 

To lowest order in Z -1 , the solutions are (t — ie)^ times Bessel functions of imaginary order 
v = i a/3/20. As for the tachyonic modes, an incoming negative frequency mode evolves to 
an outgoing negative frequency mode and there is no particle production. We solve (I6.18P 
to next order in Z _1 in Subsection 16.31 

Finally, let us consider the gauge bosons that do not acquire a mass from $i. Even 
with the deformation, there is no coupling between these light gauge fields and $i in the 
classical Lagrangian. The rolling field generates no gauge current, and therefore does not 
directly source the gauge fields. However, a coupling does arise at one loop order. As can be 
seen from the results of Ref. [58] for example, integrating out the massive Higgs, Fermi and 
gauge bosons results in an effective Lagrangian for the remaining massless gauge bosons of 
the form — \F 2 (l + a J2i=s fv ^^(M 2 /^ 2 ) j , where a = g 2 / (167r 2 ), /i is a renormalization 

scale and the trace is taken over scalars (with cs = |), fermions (with cp — 3 per Dirac 
fermion) and gauge bosons (with cy = — y). The mass matrices for the heavy particles 
depend on the orientation of the Higgs field $1 in SU(N) space. In general position, the 
simplest case to consider, the gauge group breaks to [/(l)^ -1 and all but N — 1 gauge, 
Higgs or fermi particles acquire masses of order g<p. In the undeformed theory, there are six 
scalars, two Dirac fermions and one vector so the correction vanishes by supersymmetry. 
The deformation alters the scalar mass matrix, Mj ~ g 2 <p 2 — > g 2 cf) 2 + A<^0 2 , up to numerical 
factors, leaving the other terms unchanged. The effective massless gauge boson Lagrangian 
is then — \F 2 {1 + C7 _1 ), with C a constant. The second term can also be derived from the 
finite Feynman diagram shown in Figure [TTJ 

The equation of motion for the transverse massless gauge bosons acquires a damping 
correction term as follows, to leading nontrivial order, 

A± - Cr 2 , 1 . A ± + k 2 A ± = 0. (6.19) 
(t - it) 
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Figure 11: Feynman diagram for the one-loop correction to the kinetic term of the massless 
gauge bosons. 

Ignoring the Z -2 term, the solutions are trivial and a negative frequency incoming mode 
evolves to a negative frequency outgoing mode. We solve (16.191) to order Z~ 2 in Subsec- 
tion [631 

In summary, in this subsection we have derived the equations of motion for the quadratic 
fluctuations which cannot be integrated out and which suffer significant particle production. 
In each case, to zeroth order in Z , the incoming negative frequency mode is proportional 
to {t — ie)^Hi, (k(t — ie)) and it evolves to an outgoing negative frequency mode, so there is 
no particle production across the singularity. Similar behavior was encountered in Ref . [59] , 
which studied quantum fields on compactified Milne spacetime. In that case, e was intro- 
duced to define a particular choice of continuation across the singularity. Here, we are in a 
better situation: e is a physical parameter determined by the final state of the homogeneous 
"background" and for e 7^ the time evolution of the linearized perturbations is completely 
unambiguous. Later in this paper we shall argue that, within a certain parameter regime, 
except for a small band of final states very close to the real classical solution, with e nearly 
zero, the backreaction due to particle creation is perturbatively small. 

In the next three subsections, we compute the particle production to leading nontrivial 
order in the asymptotically free coupling A^. 

6.3 Production of 5<p Excitations with a Running Coupling 

In this subsection, we shall compute the evolution of the inhomogeneous modes when the 
logarithmic running is included in an expansion in Z _1 . It is convenient to ignore e at first: 
having obtained the solution of (16.151) we can then simply replace t with t — ie and continue 
unambiguously from negative to positive times. 

The effect of the logarithm on the linearized perturbations at long wavelengths \kt\ <C 1 
and at leading order in Z -1 is easily seen. A constant time delay £ s (x) in the classical 
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background leads to a perturbation of <p of the form 



1 )~*X~lih'?%~ltt.(x)t- 2 (6.20) 
X 2 xf 

to leading order in Z -1 as t — > 0. Whereas a local perturbation of the Hamiltonian density, 
57i(x, 0) ~ p(x), leads from (14.201) to a linearized perturbation 

<50 (2) ~ p(x)0 / 0- 2 ~ /"5p( x )t 3 , (6.21) 



again to leading order in Z" 1 as t — > 0. Here, we used oc —x/x 2 — ~^ h 2 Xi/x 2 ~ Zat" 2 
at leading order in Z~ x . As we shall see, we need to include these logarithmic corrections in 
order to compute the production of particles and stress-energy perturbations as we follow 
the system across the bounce. At zeroth order in Z _1 , the particle production vanishes - the 
in vacuum maps to the out vacuum at this order. But at first order in Z _1 , there is particle 
production as well as creation of a nearly scale-free spectrum of linearized stress-energy 
perturbations. 

The physical evolution of modes on different wavelengths may be pictured as follows. 
When modes of a given wavenumber k are inside the effective horizon, i. e., \kt\ 3> 1, they are 
oscillatory and one must study their evolution using quantum field theory. However, once 
they have passed into the regime where \kt\ <C 1, spatial gradients become unimportant and 
the evolution is ultralocal. At this point, one can switch to a local, quantum mechanical 
description in which the field at every spatial point evolves independently, and in which, 
after carefully specifying the self-adjoint extension, the quantum wavefunction may be 
matched across the singularity. 

From (16.201) and (16.211) . and taking into account the corrections in inverse powers of Z, 
we expect to obtain one solution to (16.151) representing a time delay perturbation, of the 
form 

<50« = l*fV>(kt) + r? g M(kt) + . . . , (6.22) 

where p 1 ' oc t~ 2 as \kt\ — > and subsequent terms in the series are suppressed by further 
powers of Z -1 . We also expect a second solution, representing a local perturbation to the 
Hamiltonian density, of the form 

^(2) = rlfW(kt) + rV 2) (fc*) + ■ ■ ■ > ( 6 - 23 ) 

with oc t 3 at small \kt\. Substituting these ansatze into (I6.15p . /W'( 2 ) obey (16.161) . as 
in the previous subsection, which is again solved by (I6.17p . 

We can correct both these mode function solutions to next order in Z -1 by substituting 
( 16.221) and (I6.23P into (16.151) and equating the coefficients of the first subleading power in 
We find 



t 2 t 
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// 2, -t [k 2 ~l\g {2) = 3^-^, (6.24) 



t 2 y t 2 

whose particular solutions are 

9 {1) = ^W^H^ + ^in^/W^+^ + i^, 

S (2) = ^i(2H) / «(H)-icm(2H) / ( 2 )(H) + i^-^, (6.25) 

where it is conventional to define the integrals 

'1 — cosx) 



Si(y) = / dx , Cin(y) = / dx- 

Jo x Jo 



x 



(6.26) 



In fact, the solutions ( 16.251) are only determined up to an arbitrary multiple of f^ 1 ' and f^ 2 \ 
respectively (note that only one of them is consistent with the parity of either solution), 
but such an extra contribution shall not be important to us since it does not contribute to 
the particle number (n) at the leading order in l^ 1 to which we calculate. 

Equations ( 16. 22ft and (16.231) . with t replaced by t — ie throughout, solve (16. 15ft in our 
complex classical background. Since the background avoids the singularity at x = 0, there 
is no ambiguity in the evolution of the perturbation modes. We now want to take the linear 
combination of 5<p^ and Scj)^ that behaves as the incoming negative frequency mode, and 
evolve it to a time well after the quantum "bounce," when it is undergoing oscillations with 
an adiabatically changing frequency once more, and it makes sense to calculate the particle 
number. 

We work to leading order in an expansion in inverse powers of the logarithm, which we 
assume is large. The classical background solution which dominates is the "image" solution 
starting behind the singularity: as explained in Subsection l5.4[ we must define the logarithm 
to be real below the positive real x ax is (where x = 1/0) , i.e., as | (hi(x) + ln(— x) — i^)- 
Recall that at leading order in the large logarithm, i.e. ignoring \n((j)/NM) where it occurs 
inside the logarithm, we have x oc t — = t in the classical background. 

We start by extracting the time-dependence in the logarithm, as follows: 

\n((j)/NM) = I -> ln(-^r) - ^ (ln(-H) + \n(kt) - m) = l - ^ (\n(-kt) + ln(Art) - in) , 

(6.27) 

where we ignored \n(4>/NM) inside the logarithm. Our asymptotic "in" and "out" regions 
will be defined where \kt\ is large, so that the mode evolution is oscillatory, but \Mt\ is 
small so that I is still large. We now identify the negative frequency mode in an expansion 
in 

_i 

At zeroth order in Iq , the incoming negative frequency mode is just — 3/ 2 5<p^ + 

3i/ o 2 (50 ( ' 2 - ) and this evolves across t = to the negative frequency outgoing mode. At next 
order, we must include the corrections and derived in (I6.25p . When taking the 
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large \kt\ limit, we use 

n , , n 77 f°° dx sin x _ , , , , . , . f°° dx cos x 

Si(2Jfei) - - / , Cin(2fct) -> ln(fct) + 7 + / , (6.28) 

where the integrals are easily expanded as a series in l/(kt) by repeated integrations by 
parts. Notice that in the complex £-plane, Si(2&£) is an odd function of £ whereas Cin(2A;£) 
is an even function. 

For large negative times —(k/M) <C tfc < -1, we find 

-I / 1 7T 

+/ 2 I -- cos kt Qn(-2kt) + 7) - — sin fc£ 

I / 1 7T \ 1 -I 

= / 2 f --cosH - ^7-sinfctJ - -Z 2 cos H^vr + 7 + In 2). (6.29) 

Notice that the time-dependent logarithm from the first term has canceled with the loga- 
rithm in the correction to the second term. A similar cancellation occurs in 

50 (2) -> Z~* Q sin Art + 

cos fc£ I Z 2 sin kt(in + 7 + In 2). (6.30) 

Similarly, at large kt we obtain 

1 / 1 n \ 1 _x 

x) = Z 2 --cos fct + — -sinfct - -L 2 cosfct(7 + ln2), (6.31) 
\ 3 12to / 6 

and 

54> {2) -> Zq^ (-sin Art ^-cos/c£ J - -Z ^ sin /c£(7 + In 2). (6.32) 

\3 12<o / 6 

We now replace £ by £ — «e and read off the incoming negative frequency mode 

R in } 

= e" fo (-3Z~ + 3z7* <ty< 2 >) - e-% 1 + 2 ^^j (-3Z~ ^ ~ M^S^ 
-> e ifct - (k/M) < Art < -1 

-> e ifc * - ™-Me e -ikt 1 < Art < fc/M. (6.33) 
The last term describes the particle production, from (16. 8p we have 

(n) = ^e" 4 ^, (6.34) 

f n 
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Figure 12: Numerical calculation of particle production ((n) versus k) in the complex 
background of Figure 9, compared with the analytical calculation given in (16.341) . 



explicitly revealing the exponential UV cutoff at large k due to the complex classical solution 
for the background missing the singularity at \ = by the finite quantity e. 

As a check of the validity of our approximations, we have computed (n) numerically, 
using the full complex solution for the background, shown in Figure [9j and a numerical 
solution of the mode equations at each k with the field fluctuation in the incoming negative 
frequency mode. Figure [T21 compares the analytical result (I6.34p . shown as the solid curve, 
with the numerical result of the full computation. 

We can now write down the energy density in 8cj) particles produced: 

r H 3 k 7T 2 

^ = Jw k &m? ' (6 ' 35) 

where, as we have emphasized, the integral is only to be taken over those values of k so large 
that the transition from oscillatory to ultralocal evolution occurs while the homogeneous 
mode was accurately following the zero-energy scaling solution. 



6.4 Production of Light Higgs Particles 

To compute the quantum production of light Higgs particles, namely the excitations of 
$i, i = 2, ... 6, which commute with $i and whose mass is proportional to A^, we must 
solve (16.181) to first order in l" 1 . We proceed similarly to the previous subsection. Setting 
i = t — ie, we set 

$ i = l fi (f$)+g®r 1 + ...), (6.36) 

and solve (I6.18P order by order in I -1 . The two linearly independent solutions for / at small 
t are easily identified: f± ~ i^ ±ta , with a = a/3/20. At next order in (I6.18P yields 
P± = =F2i/vT5. These solutions are extended to large times as follows. Solving (16.181) to 
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lowest order in Z 1 , we find f± ~ t?J± ia (kt). At next order, we obtain 



2 , 7,2 2 A± A A± + 5 , 

9± + t^9± + k g± = -~-f± w^f±' ( 6 - 37 ) 

ot z t t A 



The right hand side is easily evaluated using standard Bessel function relations to be 
—1(5±kt~ 1 l 2 J\±ia{kt) . We now solve (16.371) using a Green's function, obtaining 



9±{t) = ±^r^& \ J±ia{hi) jf dt'J 1±ia (kt')J Tia (kt') - J Tia (kt) jf dt'J 1±ia (kt')J±ia(kt' 

(6.38) 

where we used the Wronskian t 1//2 (^j ia (kt)J^ ia (kt) — j^ ia {kt) Ji a {kt) J = (2isinha7r)/7r. The 

first term in (16. 38ft is proportional to f±{t). Using the large-argument asymptotic form for 
Bessel functions, 

/ 2 7T 7T 

J v (kt) ~ J — cos(kt-jy-- -), H>1, (6.39) 

we see the first (second) integrand in (16.381) tends to =p(sinha7r)/(7r/ct') plus oscillatory 
terms which integrate to a constant. The first integral in (I6.38P thus diverges logarithmi- 
cally, giving a contribution to g± of the form (3±hi(ki)i l l 2 J±i a {kt). This is just what is 
required to cancel the time dependence of the overall coefficient I 13 in (16.361) to first order 
in \n(kt)/lo, where Zo = ln(/c/M) and we are interested in the regime 1 <C \kt\ <S k/M as 
in the previous subsection. Recalling that Z « l — ln(— kt) + m for —k/M « M « -1 
and Z « Zq — ln(A;t) for 1 <C Zct <C fc/M, we can check that the same cancellation occurs 
for large negative t using J_ ia (/ct) = e na J-i a (—kt) and J\ + i a {kt) = — e~ na J\+i a (— kt), which 
hold when t is analytically continued to negative values by passing below the origin. The 
coefficient of J± ia (kt) in the first term of (16.381) is hence an even function of t, tending to 
/3±ln(|Atf|) for large \kt\. 

The second term in (I6.38P involves an integral which converges at large \t'\. Here, we 
can make use of the formula / °° dxJ\ + ^{kx) J^(kx) = l/(2k), [60] to obtain for our two 
linearly independent modes, to first order in Zq 1 and at large \t\ (still <C M _1 ), 

$i,± - le^l^i-t) 1 ' 2 (j±ia(-kt)(l + ^) =F . m ^ J^U-ki)) , kt « -1 

\ Zo Z sinn arc J 



-> lot 1 J±ia(kt) T , J Tia (kt) , kt > 1. (6.40) 

\ Z smh a7r / 

Using (I6.39P we now construct the incoming negative frequency mode as a linear combi- 
nation of $i,±, and then use (I6.40p to evolve it to large positive kt. Explicitly, at large 
negative kt, we have 



R+ ~ e ifet ~ e-^e-^J f J fa (-fct)(- - -) + J- ia (-kt)(- + -)) , « -1, 

2 V 2 \ cs c s J 

(6.41) 
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where c = cosh(7ra/2) and s = sinh(7ra/2). Evolving to large positive kt using (16.401) and 
working to first order in Iq 1 , we find 

R+ ~ (e lkt (l + 2i^-) + % e~^e~ 2 ^l± cothvra^ . (6.42) 

\ *0 H) / 



The last term yields our final result for the mode-mixing coefficient. Using (3 + = — 2i/vl5, 
a = a/3/20, we infer the energy density in light Higgs particles 

J (2tt) 3 151n 2 (fc/M) ' V ; 

which is parametrically of the same order as the result (I6.35f) for 8(f> particles. However, 
we have still to multiply by the total number of light Higgs particles. The simplest case to 
consider is where the SU(N) symmetry is broken by $1 to [/(l)^ -1 . In this case, there are 
of order N surviving light Higgs bosons whose mass will be due solely to the deformation 
term. Hence the result (16.431) is multiplied by N. (In more general cases, there could be of 
order N 2 such light Higgs modes.) 

6.5 Production of Massless Gauge Bosons 

To calculate the production of remaining light gauge bosons, we need to solve equation 
(16.191) to the first nontrivial order in The effect of the l~ 2 correction is more modest in 
this case and a suitable ansatz for the transverse mode functions is 

A ± = A sin kt (1 - Cr 1 ) + A 1 cos kt + gl~ 2 + .... (6.44) 

Substituting into (16.191) we obtain 

g + k 2 g = AoC^kr 1 cos kt - t~ 2 sin kt) + A x Ckt~ x sin kt, (6.45) 

whose solution is 

g = \ CA v [ sin kt (Cin(2H) - 1) + cos kt Sin(2fct)] 

+ \CA 1 [sin kt Sin(2kt) - cos kt Cin(2Jfet)] . (6.46) 



Using the large time limits of the Cin and Sin functions, as given in Subsection I6.3[ we may 
obtain the behavior of an incoming negative frequency mode A^' m ~ e lkt , k/M <C kt ^ 
— 1, including the effect arising from the continuation of /, to obtain 



A (+),in _^ e m + ™C e -*t. ( 6 .4 7 ) 
2l 



From this we read off the mode mixing coefficient, and hence the energy density in created 
light gauge bosons: 

^S§f k ^mr'-' (6 ' 48) 
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which is suppressed relative to the other cases by the square of the logarithm. 

As for the light Higgs particles, we can discuss the dependence on N. If the SU(N) 
symmetry breaks to [/(l)^ -1 , then there are of order N massless bosons remaining, so the 
result (I6.48P is multiplied by N. However, the diagram which gives rise to the coupling C 
involves only one closed line (hence one power of N) and is proportional to A^ ~ iV~ 2 Z -1 . 
Hence C scales as N . (For symmetry breaking to a larger group, there can be of order 
N 2 massless gauge bosons remaining.) The quantum production of massless gauge bosons 
is always subdominant to light Higgs production, by two powers of l/ln(fe/M) and by two 
powers of 1/N. 



6.6 Backreaction 

We now want to estimate the backreaction effect of particle production on the homogeneous 
mode of the field <f>. We do so using energy conservation: the total energy is always zero 
so any energy which goes into particles will lower the energy available in the homogeneous 
mode and therefore cause it to sink lower down the negative potential. 

One possible criterion for the backreaction to be small is that the energy density p c in 
created particles should be small compared to (minus) the potential energy density reached 
at times \t\ ~ e. In the previous subsections, we have shown that most of the energy density 
goes into light Higgs particles. Keeping in mind that the number of light Higgs species is 
of order N and that for most of the final wavefunction e ~ W 3 RAds/(N\ ki(MRAds)\ )> 
with a minimal spread wavepacket corresponding to W = 1, we find from (16.431) that 

N N 5 
Pc ~ e±\HMR AdS W ~ W^R\ dS (6 ' 49) 

Comparing this with (minus) the potential energy at \t\ ~ e, 



we see that the energy density in created particles is smaller by powers of 1/N and 
l/\\ia(MRAds)\, suggesting that backreaction is small. (We are assuming MRAds <S 1, 
so that the quartic coupling near the starting value of <fi, f ~ 1/| hx(MRAds)\i is small and 
perturbation theory reliable.) 

A more conservative criterion for the backreaction to be small is that the homogeneous 
mode should keep enough energy to roll back up the hill to nearly its starting value. The 
potential for is 

Setting V = to obtain the starting value of 0, we find 
T 1 



' start 



2,X Q R A dsH<Pstart/NM)^ V 3 A RAdS I ln(MR Ad s) 1 1/2 ' 



(6.52) 
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The minimal value of (f> reached after the bounce is obtained by setting V((f> m i n ) = —p c . 
Demanding that <p m i n — <fi s tart *C <fi s tart i n order to have small backreaction, one finds the 
condition p c <C N 2 \ hi(MRAds)\l R\ds'i usrn g (16.491) . this becomes 



backreaction is small as long as | bi(MRAds) I ^ 1) f° r final values of <ft carrying most of the 
probability. This value of W corresponds to a spread Axi ~ iV~ 3 / 4 | hi(MRAds)\~ l ^ 2 RAdSi 
parametrically small compared to the distance RAdS traveled in %-space during the bounce. 
Hence our wavepacket is well-described by a classical bouncing solution, with small backre- 
action. The fraction of the probability associated with final values for which e is so small 
that backreaction is large (and hence our calculation breaks down) is ~ | ln(Mi? J 4ds)| _1//4 • 
Hence we conclude that when the logarithm is large, the most probable outcome is a bounce. 

In Appendix |Bl we mention that the beta function for the coupling / of our double 
trace deformation is of order f 2 and does not receive higher loop corrections in the large N 
limit. However, in the large N limit we considered in the appendix, the 't Hooft coupling 
g t as well as / were kept fixed. In the present situation, if we choose a minimal width 
wavepacket, for which we have to choose / <C 1/iV 3 in order to have small backreaction, we 
have to worry about subleading terms in the large N expansion. For instance, if the beta 
function contained terms linear in /, they might interfere with the logarithmic running of 
/ that led to its asymptotic freedom. In Appendix (Bj we show that the one-loop effective 
potential does not receive contributions proportional to fg t . It is conceivable that higher- 
loop diagrams do give contributions to the beta function for /, for instance of order fg 2 /N 2 . 
If / <C 1/iV 3 and if g t were kept fixed for large N, such contributions would dominate the 
f 2 contribution considered in this paper. Therefore, for a minimal width wavepacket we 
would have to choose gt smaller than the appropriate power of 1/N, in which case our 
analysis would go through. This would deepen the gap between the regime in which our 
computations are reliable and the large 't Hooft coupling regime in which the bulk theory 
can be reliably described in supergravity. Either way, it will be interesting to explore in 
future work to what extent our results can be extended to large g t . In any case, it appears 
that such problems can be avoided by choosing a sufficiently broad initial wavepacket, as 
described in the previous paragraph. 

Now we would like to tie up some loose ends from previous (sub) sections. First, are 
non-linear corrections to the fluctuation equation (16.151) really small? It is easy to see from 
the full equation of motion that non-linear corrections will be small as long as \5<j)\ <C \<p\. 
In our setup, the perturbations 5<p have a quantum mechanical origin: as the singularity 
is approached, more and more Fourier modes "freeze" and start behaving classically. To 
estimate the typical magnitude of the perturbation 5(f>, let us compute the contribution 
of frozen modes to the expectation value of 5(f) 2 . From (16.171) . we see that a properly 




(6.53) 



So, for example, if we choose 



W ~ iV 1/4 , 



(6.54) 
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normalized frozen negative frequency mode with wavenumber k behaves like 

(6.55) 



gik-x 



so that 



1 f d 6 k H 2 

(50(x)50(x))| frozcn ~-J ( Ll~ f™s } (6.56) 

with t = t — ie (see Subsection 6.2) and where in the last step we have used the infrared 
cutoff k > 1/Racis due to the finite volume of S 3 . So the typical size of fluctuations at time 
t is 



12 

This we need to compare to the homogeneous background 

1 



(6.57) 
(6.58) 



The most dangerous region is around t = 0, where 

|t| = e ~ W 3 y /XiR AdS . (6.59) 

Therefore, nonlinearities will be unimportant as long as W> 1, which is satisfied for our 
choice ( 16341) . 

Second, will the fluctuations have a significant effect on the homogeneous mode? In 
other words, was it a good approximation to first solve for the homogeneous mode, setting 
the fluctations to zero, and then solve for the fluctuations in this homogeneous background? 
In this Subsection, we have already performed a consistency check on this, by verifying that 
the homogeneous mode keeps enough energy to bounce back almost to its starting value. 
However, one could ask about effects at intermediate times, in particular near t = 0. 
Again it is easy to see from the full equation of motion that effects of fluctuations on the 
homogeneous mode will be small as long as \5<p\ <C \<fi\. So the analysis of the previous 
paragraph applies and our approximations are good as long as W ^> 1. 

Finally, in Subsection 4.2 we assumed that the surface £ where the classical field 4> 
reaches infinity, given by t = t s (pc), was smooth and space-like. Now we would like to verify 
that |Vt s | < 1 so that E is indeed space-like, for typical classical solutions <fi generated by 
the quantum mechanical initial state we have considered. Below (14.191) . we concluded that 
t s is the time delay mode. In the linearized approximation, close to the singularity and to 
leading order in gradients, it can be related to the field fluctuation 8(f): 

t s (x) = -t\t=o ~ V^P<ty|t=o- (6.60) 



To estimate Vi s , we compute the contributions of the frozen modeo to the two point 
function 

/k 2 f rlk 

d3k ^i ~ X ^J T- (6 ' 61) 



14 These are the only modes that have become classical and thus contribute to the classical perturbation 
that determines the shape of E. 
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This is indeed small compared to 1 since A^ is very small and the integral has both an IR 
cutoff due to the finite volume of S 3 , and a UV cutoff because modes with k > 1/e never 
freeze: we conclude (16.611) is suppressed by inverse powers of the large parameters N and 
\a(MRj\ds)- Thus typical surfaces £ generated from our quantum initial conditions are 
indeed smooth and spacelike, as assumed in Subsection 4.2. 



7 Stress-Energy Correlators 

In addition to computing the quantum creation of particles, we would like to calculate the 
two-point correlation function of the stress-energy tensor for the scalar field fluctuations as a 
function of time. There is no gravity (and hence no backreaction) in our system, but stress- 
energy tensor correlators of the boundary CFT determine the metric perturbations in the 
bulk, i.e., the cosmological spacetime, and are as such related to cosmological perturbations. 
In fact, the relation between the stress-energy tensor and the bulk metric is completely 
analogous to the relation between the operator O and the bulk scalar field in Section [31 
We shall consider the conformally improved stress tensor [61], consisting of the canonical 
stress tensor plus an identically conserved quantity which has been added to make the stress 
tensor traceless in the limit when the theory is conformally invariant. In our case, this limit 
is where the logarithm is large so the running of A may be neglected. In four spacetime 
dimensions, the conformally improved stress tensor is 

T, v = (d^d^ - g,u{\{d<P) 2 - ^A)J - l - {V,d v - g, u d 2 ) 2 . (7.62) 

We shall calculate the correlators of T^ v to linear order in the field fluctuation 50 about 
the background solution. We shall treat the latter as a classical background, although 
strictly speaking, as we detailed in Subsection l5.6[ the expectation value of the homogeneous 
background mode does not exist. Nevertheless, at times well after the bounce, in the 
parameter regime where backreaction is small, we expect the description of in terms of 
a localized wavepacket rolling back up the hill is accurate, so we can replace with the 
appropriate classical solution and take expectation values of the inhomogeneous fluctuation 
modes £0(x) using the Schrodinger wavefunction, calculated above. As we shall discuss 
below, these linearized correlators are in principle sufficient to determine the correlators of 
the linearized metric perturbations in the bulk. 

Denoting the linearized perturbation of the stress tensor as 5T flu , the linearized compo- 
nents at given wavenumber k are: 



• 2 d ( 5<j)\ 1 2l 



5T 0i = iki ( -050 



3 r T 3 

5T tJ = 5 l3 ( ^050 + jj<M0 - ^050 - ^<j>k 2 6<t> ) + ^0^50. (7.63) 
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We want to compute the fluctuations in the stress tensor in the incoming adiabatic vacuum 
state for the fluctuations, denoted |0, in), relative to the vacuum defined by the outgoing 
adiabatic vacuum state |0, out). Equivalently, we compute the correlators of T^ u normal 
ordered with respect to |0, out). To compute the correlators, we first express the scalar field, 
and the linearized stress tensor, in terms of the adiabatic positive and negative frequency 
modes, and creation and annihilation operators, appropriate to |0, out). We then normal 
order the expression and calculate the expectation value in |0, in). This requires expressing 
the "out" creation and destruction operators in terms of "in" creation and destruction 
operators: 

a out = aa in + (3*a[ n , oJ ut = a*a\ n + (3a- m . (7.64) 

Consider some quantity linear in the quantum field fluctuation and its time derivatives: it 
may be expressed as 



W> X ) = E "TtW (^txittV^ + al outX ( ~le~^) , (7.65) 
where xiut"' 1 are the positive and negative frequency parts with respect to the out vacuum. 



We normalize the creation and destruction operators so 
vacuum-subtracted expectation of some operator O as 



Ok, a k , 



5 k k ' etc. We define the 



(O) = (0,in|C|0,in) - (0, out|O|0, out). (7.66) 

Substituting ( 17.641) into ( 17.651) and ( 17.66[) . using the commutation relations for and a k 
and their action on the "in" and "out" vacua, and replacing ^ k with V J d 3 k/(2n) 3 , we 
obtain the two-point correlator 

(/(f,x)/(f,0)) = J (l^^^^ (2|/?| 2 |X^| 2 + ^^(xi^) 2 + «/3*(x^ } ) 2 ) , (7.67) 

where we used \a\ 2 — \f3 2 \ = 1. In the regime of interest, where f3 is small, we can replace 
a with —1. 

Since we only want the leading order result, we can use the positive frequency mode 
functions to zeroth order, namely 



v kt (kt) 2 y 

which leads to the following positive frequency contributions to ST^: 

^<t> rrp(+) _ k 2 -ikt 

2_ 6I ° - 3t 6 

X A ST M _ ik 1 + M -ikt 
2_ 01 ~ 1 3t 2 
% (+) _ l + ikt 1 / 3i 3 \ —ikt (7Qg) 
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Notice that the spatial trace 8T^' = 8T^ as it has to be from the tracelessness of the 

stress tensor. Hence it is only the traceless part ST^' = ST^ — j-SijST^ which needs to 
be considered: 




2 H% = [t,k, - -k%) - [1 - E - w j (7.70) 

Equations (17.691) and (I7.67P determine the two point correlation function (at equal times) 
of the stress-tensor in the boundary theory. Here we shall present only a partial discussion 
of the equal time correlators, postponing a fuller discussion to future work in which we 
attempt to use the boundary correlations to determine the fluctuations in the bulk. 

The first term in (17.671) provides the simplest contribution. In order to compute it, 
our strategy is to replace factors of k 2 under the integral by minus the spatial Lapla- 
cian outside the integral. We then perform the angular integrals in (17.671) to obtain a 
factor of 4n(smkr)/kr. The two powers of k in the denominator then cancel the k 2 
factor in the measure, leaving us with r _1 J °° dksmkr. When performing the integral 
over \/3\ 2 , we replace k with r _1 , the value of k dominating the integral. For example, 
J °° dk(sinkr)/ (ln(fc/M)) 2 « r -1 / (ln(l/Mr)) 2 . These manipulations are sufficient to com- 
pute the first term in (I7.67p . The second and third terms involve additional oscillatory 
factors of e ±2lkt under the integral, which generates additional structure on the causal scale 
r = 2t, but which suppresses them for r<t. 

In this way, we obtain the equal time correlators of various scalar quantities constructed 
from the stress-tensor, for example we exhibit the detailed form for r Ci, 

\n(l/Mt) 2 



\ 7i 2 (\nMr) 2 3rH 2 

ln(l/Mt) / 20 2 \ 



(5T 00 (r,t)5T w (0,t)) 

{di 5Mr,msT 0i (o,t)) = A ;;— )2 {^-^) 
m^ dkdM - J» (i- + + ( , 71) 

The key point is that all of these correlators are nearly scale-invariant in form, that is 
they are all of the form t~ n g(r/t), where the overall power of t is determined on dimensional 
grounds. The only deviation from perfect scale-invariance is due to the logarithmic running 
of the coupling with scale. As discussed in the introduction, this behavior is the result of 
the classical scale-invariance, and the modest quantum breaking of scale-invariance, in the 
boundary conformal field theory. 

Although the background energy density is zero, we can define a dimensionless pertur- 
bation by comparing the fluctuation 5Tqq = 5p with the background quantity p + P = <fi 2 . 
We find for example 

5p(r,t) 5p(0,t) 1 f_ 

{ P + p P + p ' ~ iV 2 ln 2 (l/Mr)ln(l/Mt)r 6 ' K ' ' 

a slightly red spectrum, with the reddening arising from the running coupling as before. 
Furthermore, we note that the overall amplitude of the fractional energy density pertur- 
bation is proportional to l/l 3 , where l/l is the appropriate value of the running coupling, 
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r=0 boundary 



Figure 13: The stress-energy correlators on the boundary may predict the spectrum of 
classical cosmological fluctuations on a spacelike surface in the bulk at late times. 

which has been our small parameter all along. To summarize, the fluctuations in give rise 
to stress-energy perturbations that are naturally small, being suppressed by powers of 1/N 
and l/l, approximately scale-invariant, with a slightly red tilt due to asymptotic freedom, 
scalar, adiabatic and nearly Gaussian in character, to leading order in l/l and 1/N. 

The quantities of physical interest, however, are correlators of metric fluctuations in 
the bulk cosmology, which couple to the stress-energy tensor of the boundary theory. We 
would like to compute the bulk perturbations on a spacelike surface at late times, roughly 
halfway between the past and future singularities in our model, as shown schematically 
in Figure [131 In general, while knowing the state of the boundary theory is in principle 
equivalent to knowing the state in the bulk theory, it is in practice not straightforward to 
extract local information about the bulk theory from the boundary theory (see [62] for a 
recent discussion). Our specific problem might be simpler due to the fact that the metric 
fluctuations behave classically in the relevant regime. Furthermore, it is sufficient to treat 
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them in a linearized approximation. 

One idea for how one might be able to compute the spectrum of bulk metric perturba- 
tions is as follows. The state of the boundary theory is encoded in the boundary correlators, 
some of which we have computed in this paper. Since the relevant modes are frozen in, 
it may be possible to view the boundary correlators as those of a classical ensemble of 
fluctuating fields. Each member of the ensemble could then be viewed as an expectation 
value of a boundary operator, which would determine a corresponding bulk solution (see for 
instance [46]). For instance, knowing the expectation value of the stress-energy tensor on 
the boundary determines a solution for the bulk metric. So we would end up with a classical 
ensemble of bulk metric perturbations, which we could use to compute bulk correlators. 

Furthermore, since we are interested in wavenumbers k 3> R^lsi which is the only scale 
that enters in the correspondence, and since scale-invariance is holographic in the sense 
that the fluctuations on a planar subspace of a space with scale- invariant fluctuations are 
scale-invariant, it seems plausible that the bulk metric perturbations will turn out to be 
scale-invariant too. 

8 Conclusions 

We have studied the AdS/CFT dual description of Af = 8 gauged supergravity in five 
dimensions where smooth asymptotically AdS initial data can evolve to a big crunch singu- 
larity in the future. At the singularity the classical supergravity description breaks down, 
but we have argued that the dual field theory evolution provides a consistent quantum 
description of the singularity. Specifically the boundary theory predicts a transition from 
the big crunch to a big bang as the most probable outcome of cosmological evolution. 

If our results extend to more realistic cosmologies this would place cosmological models 
involving a big crunch/big bang transition on much firmer footing. It is tempting indeed 
to view AdS as a maximally symmetric box describing a smooth, isotropic patch of the 
universe which, in the case at hand, collapses to form a spacelike singularity. This applies 
for example to the cyclic universe models where the universe collapses after a phase of dark 
energy domination rendering it locally Minkowski on scales corresponding to a Hubble 
volume today. Since in the near singularity regime the cosmological constant appears to 
play no significant role, one would then use AdS/CFT, with boundary conditions that are 
a minimal deformation of the maximally symmetric setup, to model a patch of the universe 
like this through the singularity. The methods we have developed here may therefore allow 
the AdS / CFT correspondence to be used as a "laboratory" in which one can construct and 
study a large variety of cosmological scenarios. 

An important limitation of our work is that we have good control over the field theory 
calculations only at small 't Hooft coupling, for which the bulk is in a stringy regime. There 
are indications, however, that our results may extend to strong coupling. In particular, the 
quantum effective potential for the unstable scalar in the boundary theory is negative and 
unbounded also at large 't Hooft coupling. Therefore the fact that AdS/CFT maps the 
problem of bulk cosmological singularities into a field theory instability of this type is a 
universal feature of this class of models. 
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We should also note there is a narrow band of the final wavefunction in the boundary- 
theory, centered around the value predicted by the real classical bouncing solution, that 
cannot be calculated in perturbation theory. For final values within this band, which 
represents a parametrically small fraction of the total probability, we do not know yet 
whether the dual theory predicts a bounce. 

Finally it would be interesting to find a mathematical framework, describing unstable 
field theories like the one we have been discussing, where the arbitrary phase associated 
with the choice of self-adjoint extension is fully determined. A priori the phase leads to an 
arbitrary function when one generalizes the self-adjoint extension to quantum mechanics in 
higher dimensions and to field theory. The self-adjoint extension, however, should really be 
viewed as part of the definition of the theory - remember that without a boundary condition 
at large field values, one simply cannot discuss the quantum evolution of a wave packet - so 
it is natural for it to have the same symmetries as the Lagrangian. This, together with the 
ultralocality of the evolution, reduces the arbitrariness again to the single free parameter 
present in one-dimensional quantum mechanics. 

A remarkable consequence of the instability and approximate scale-invariance of the 
boundary theory is that an approximately scale-invariant spectrum of stress-energy pertur- 
bations is generated on the boundary, whose amplitude is suppressed by the asymptotically 
free coupling that governs the instability. If it can be shown that this corresponds to nearly 
scale-invariant bulk metric perturbations, and if our setup can be generalized to more re- 
alistic models, this might provide an attractive mechanism for explaining for the origin of 
the observed density variations in the universe, since the approximate scale-invariance of 
the perturbations arises directly from the underlying approximate conformal invariance of 
the dual field theory. 
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A More on the Bulk Theory 



In this appendix, we first show that the bulk theory discussed in Section 2, with action 



l -R - i(V0) 2 + — L- (I5e 2 ^ + 10e" 4 ^ - e" 10 ^) 

2 2 4i MdS 



admits an 0(5)-invariant Euclidean instanton solution of the form 



(A.l) 



ds 2 = R\ dS (^y + P 2 ^ 4 ) (A.2) 

with = 0(p) and for scalar boundary condition atf = f/3. 

The field equations determine b in terms of 0. Asymptotically one finds 



, i t 2a 2 (hip) 2 a(4p-a)\np 8/? 2 - 4a/? + a 1 
~V + V + 12? 



& = P 2 + 1 + 3/ + o„ 2 + , < A.3) 

and the scalar field itself obeys 



b 2 4>" + (y + bb ') <P' ~ RmsV* = 0, (A.4) 

where prime denotes d p . 

Regularity at the origin requires <f)'{0) = 0' o = 0. Thus the instanton solutions can be 
labeled by (f>o, the value of at the origin. For each O , one can integrate ( 1A.4I) and get 
an instanton. Asymptotically one finds 0(p) = alnp/p 2 + f3/p 2 , where a and (3 are now 
constants. Hence for each 0o one obtains a point in the (a,/3) plane. Repeating for all 
0o yields a curve oti{(3) where the subscript indicates this is associated with instantons. 
This curve is plotted in Figure O (right panel) . The left panel of Figure dH shows 0o as a 
function of (3. 

The slice through the instanton obtained by restricting to the equator of the 5* 4 defines 
time symmetric initial data for a Lorentzian solution. The Euclidean radial distance p sim- 
ply becomes the radial distance r on the initial data slice. So given a choice of boundary 
condition a(/3), one can obtain suitable initial data by first selecting the instanton corre- 
sponding to a point where the curve oti((3) intersects a(/3), and then taking a slice through 
this instanton. 

One immediately sees that all «/ boundary conditions, for / > 0, admit precisely one 
instanton solution. When / — ► one has 0o — > oo. The dual field theory is supersymmetric 
in this limit, so one expects there to be no regular instanton solutions that describe the 
decay of the AdS state in the bulk theory. For small positive / one sees that (3 ~ 1//. This 
scaling is nicely reproduced by the solutions of the simplified model field theory discussed 
in the text. From the asymptotic form flA.3j) of b it follows that the conserved mass (cf. 



(12.121) ) of the resulting initial data is M = -ir 2 R 2 AdS a 2 /4, so M w -2R 2 AdS for small /. 
The instantons therefore specify initial data with negative mass in this theory. It is only 
for AdS-invariant boundary conditions discussed below that instantons of this type define 
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Figure 14: The value of at the origin as a function of (3 (left), and the function ai{(5) 
(right) obtained by integrating the field equations outward starting with different values 
for O . 



initial data of exactly zero mass, in line with their interpretation as the solution AdS^ 
decays into0 

Lorentzian Evolution 

Analytic continuation of the instanton geometry yields a Lorentzian solution that de- 
scribes the evolution of instanton initial data under AdS-invariant boundary conditions 
[10,42]. These differ slightly from otf boundary conditions and are most conveniently ex- 
pressed as [38,39] 

Oi{l- f - In a) = fP, (A.5) 

where / is again an arbitrary constant. One sees in particular that rescaling r leaves / 
unchanged. One also sees that the otf boundary conditions considered in the text are 
approximately AdS-invariantor when / is small, which is the parameter regime we will 
concentrate on. This restricts us to positive / since / < boundary conditions admit 
instanton solutions only for |/| > 0(1). 

After analytic continuation the origin of the Euclidean instanton becomes the lightcone 
of the Lorentzian solution. Inside the lightcone, the 5*0(4, 1) symmetry ensures that the 
solution evolves like an open FRW universe, 

ds 2 = -dt 2 + a 2 {t)da A} (A.6) 

where da^ is the metric on the four- dimensional unit hyperboloid. Under evolution rolls 
down the negative potential. This causes the scale factor a(t) to vanish in finite time, 
producing a big crunch singularity. Outside the lightcone, the solution is given by (1A.2j) 
with df2 4 replaced by four-dimensional de Sitter space. The scalar field remains bounded 
in this region. On the light cone we have = 0o an d d t <p = (since P = at the origin 
in the instanton). 

15 One can also show that these instantons have finite action, which is large for small / and small for 
large /. 
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One can verify that this analytic continuation indeed satisfies the boundary conditions 
( 1A.5|) by doing a coordinate transformation in the asymptotic region outside the light cone. 
The relation between the usual static coordinates (12. 7p for AdS?, and the SO (4, 1) invariant 
coordinates, 

ds 2 = R 2 AdS (j^Tf + P 2 (~ dr2 + cosh2 Tdn ^) > ( A - 7 ) 

is 

p 2 = r 2 cos 2 1 - sin 2 1. (A.8) 
Hence the asymptotic behavior of in global coordinates is given by 

, . . a In r 8 , . „, 

<K r ) = — + £rn (A.9) 

where a = a^/ cos 2 t and /? is given by (IA.5I) with a replaced by a. Hence a is now time 
dependent and blows up as t — > 7r/2, when the singularity hits the boundary. 

Although for a/ boundary conditions one can in general not obtain the solution by 
analytic continuation of the instanton, we have argued in Section 2 that when / is small and 
positive, the effect of the change in the boundary conditions on the evolution is negligible 
except in the corners of the conformal diagram where the singularity hits the boundary 
at infinity. In particular, since the evolution of the initial data has trapped surfaces, a 
singularity will still form in the central region. However, a priori it is now possibe that the 
singularity is enclosed inside a large black hole and does not extend out to infinity. We now 
investigate this possibility. 

Black Holes with Scalar Hair 

We now numerically integrate the field equations derived for static spherical solutions 
to verify if the theory flA.lj) for aj boundary conditions admits a class of static, spherical 



black holes with scalar hair outside the horizon. In light of the discussion above, we would 
like to know whether there are negative mass black holes with scalar hair@ Writing the 
metric as 

dsl = R\ dS {-h(r)e- 25{r) dt 2 + h-\r)dr 2 + r 2 rffi 2 ) , (A.10) 
the Einstein equations read 



/3ft 



r 



h( f>,rr + 3^ + K ) 0,r = &AdS V fr ( A ' U ) 



r 2 2 



2(1 - h )- rh r - -^ r h = -r*R A AdS V{4>), (A.12) 
8 = - V . (A.13) 
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Negative mass hairy black holes were found in closely related theories in four dimensions [11]. 
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Figure 15: The function abh{P) obtained by integrating the field equations outward starting 
with different values for e at the horizon R e = 0.1 of spherical hairy black holes. 

We integrate the field equations outward from the horizon. Regularity at the event horizon 
r = R e requires 

*'«■>- 5rSfes- <A14) 

The scalar field asymptotically behaves as (12. 8ft . so we obtain a point in the (a, 0) plane 
for each combination (R e ,<f) e ). Repeating for all e gives a curve a.R e (/3). We show this 
curve in Fig [15] for R e = .1. Given a choice of boundary conditions ct/(/3), the allowed 
black hole solutions are simply given by the points where the black hole curves intersect 
the boundary condition curve: abh(R e , P) = a f(P)- One sees that for all values of / > 
there is precisely one hairy black hole solution of size R e = 0.1, and the same is true for 
different radii R e . For af boundary conditions the mass of the hairy black holes is given by 

M h = Q[d t ] = 2n 2 R 2 AdS 0M O + P 2 (l - i/) ) , (A.15) 

where M is the coefficient of the 0(l/r 6 ) term in the asymptotic expansion of the g rr 
component of the metric. 

This is shown in Figure [TB| as a function of horizon size R e , for two different choices of 
af boundary conditions (/ = 0.05 (top), and / = 0.1 (middle)) and for Schwarzschild-AdS. 
One sees that the scalar hair increases the mass; the hairy black holes are more massive 
than a Schwarzschild-AdS black holes of the same size. In particular, one sees that Mh > 
for all R e . For large R e one has E h ~ R 2 Ad s R i with M h/M s > 1 for all R e and M h /M s -> 1 
for large R e . Furthermore, the mass (for given R e ) increases for decreasing /. 

The dual interpretation of the increase in mass relative to vacuum black holes is that 
whereas Schwarzschild-AdS black holes correspond to thermal states in the usual vacuum 
of the dual field theory, hairy black holes are described by typical excitations around the 
local maximum of the effective potential of the dual theory in the presence of the negative 
double-trace deformation [11] (see Fig [2] (left)). An important consequence is that the 
singularity that develops from initial data defined by the instanton cannot be hidden behind 
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Figure 16: The mass M of static spherical hairy black holes (in units of R\ds) as a function 
of horizon size R (in units of Rms)i f° r two different choices of atf boundary conditions and 
for Schwarzschild-AdS. From top to bottom, the curves correspond to / = 0.05, / = 0.1 
and Schwarzschild-AdS. One sees that all hairy black holes have positive mass. 



an event horizon, since there is no candidate black hole with the required mass. Instead 
the singularity probably extends all the way to the boundary, cutting of all space. 

By contrast there do exist negative mass black holes for bulk boundary conditions that 
correspond to a deformation of the dual field theory that yields an effective potential with 
a global negative minimum, such as shown in Figure [2] (right). These negative mass hairy 
black holes have a natural interpretation in the dual theory as typical excitations above the 
global vacuum state [11]. 

For example, in the case at hand, for boundary conditions a/ |£ = f (3 — e/3 3 with e 
sufficiently small, the boundary condition curve o;/ )e will generally have two intersection 
points with the black hole curves abh{Re) of Figure [151 The first set of intersection points 
specifies a branch of hairy black holes with rather small (3 and hence with not much hair. 
Their mass is positive and they correspond to excitations around the local maximum of the 
effective potential in the dual theory, as before. On the other hand the solutions with the 
larger (3, which are associated with the second intersection point, have much smaller mass 
and provided e is sufficiently small a subset of these have negative mass. For example, for 
/ = 1 and e = .05 the theory admits an M ~ — 3R\ dS black hole of horizon size R e = 0.1 
and with cb e = 1.517. 

The existence of an additional branch of hairy black holes - some with negative mass - 
for af >e boundary conditions also means that in this case there does exist a hairy black hole 
with the same mass as the instanton, M « — eR^ds- This is therefore the natural end state 
of evolution of the initial data defined by the instanton, under the modified boundary 
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conditions. In general, regularization of the dual theory that makes the effective potential 
bounded from below encloses the singularity in the bulk inside a horizon. Furthermore if 
one considers a series of dual theories where for decreasing values of e one finds the horizon 
size of the black hole of the same mass as the instanton increases. In the limit e — > the 
black hole becomes infinitely large and we recover the cosmological solutions we studied 
here. 

B Renormalization of the Boundary Theory 

In this appendix, we discuss M = 4 super- Yang-Mills theory with gauge group SU(N), 

S = J d A xTi j-ii^F^ - ^D^D^ + ^g 2 [&, &} $,] + fermions j , (B.l) 
deformed by the double trace potential 

W = -£ ! d^xO 2 , (B.2) 



2 

see (13 .3p . In our conventions, 

Ffu> = d^A u - d v A^ + ig[A^ A v ], = d^ 1 + ig[A„, (B.3) 

The operator O was defined in (13.11) : 



O = — Tr 

N 



5 

i=2 



(B.4) 



As explained in [13] and reviewed in Section 3, the computation of amplitudes at order 
f 2 involves matrix elements of 

^ j d 4 xd 4 y0 2 (x)0 2 (y). (B.5) 

From conformal invariance, (O(x)O(y)) = v/\x — y\ A on flat H 4 . This leads to a short 
distance divergence that renormalizes / and survives in the large iV limit (where the 't Hooft 
coupling g t = g 2 N as well as / are kept fixed): 



J d 4 xd 4 yO(x)0(y)(0(x)0(y)) ~ n 2 f 2 \nA J d 4 x0 2 (x), (B.6) 



with A an ultraviolet cutoff. This leads to a one- loop beta function for /, which does not 
receive higher loop corrections in the large N limit [13]. 

Diagrammatically, we represent the (Tr$ 2 ) 2 vertex in double line notation as in Fig- 
ure [TTl The diagram in Figure [18] then leads to the renormalization discussed above. To 
see that it survives the large N limit, note that both vertices come with a factor f /N 2 , 
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Figure 17: The (Tr$ 2 ) 2 vertex in double line notation. 




Figure 18: The one-loop diagram that renormalizes the coupling / at large N. 




Figure 19: A diagram that does not survive the large N limit. 
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Figure 20: Higher order diagrams that survives the large N limit are factorizable. 







Figure 21: One- loop diagrams and a counterterm for a single scalar field. 



while the two index loops each contribute a factor of N, so that the one-loop diagram in 
Figure [18] is of the same order in N as the vertex Figure [T7] (if / if kept fixed in the large 
N limit). This also makes it clear that diagrams such as Figure [T9l do not survive the large 
TV limit: there are no index loops that can provide factors of N. By the same token, the 
higher-loop diagrams that survive the large N limit are those like the one in Figure [20l 
which are factorizable. 

To verify the above one-loop result and to show the absence of higher-loop corrections 
to the beta function for /, we wish to compute the diagrams that survive the large N 
limit. Before doing so, let us briefly review the one- loop renormalization of the theory of 
a single (massive or massless) scalar field with a — A0 /4 potential. We use renormalized 
perturbation theory (see, for instance, [63]): 



(B.7) 



where all the fields and parameters are renormalized ones. To one-loop order, the two to 
two scattering amplitude gets contributions from the diagrams in Figure [2IJ The result is 



iM(pip 2 -> p 3 Pi) = 6iA + (6iA) 2 \iV{s) + iV(t) + iV{u)} + 6i5 



(i) 

A > 



(B.8) 



where s, t and u are Mandelstam variables, 

d 4 k 



iV(p 2 



(27r) 4 fc 2 + m 2 (k + p) 2 + m 2 



(B.9) 



and 5^ is the order A contribution to 5\. We impose the renormalization condition 



iM(pip 2 -> p 3 Pi) = 6iA at {p x + p 2 f = (pi + p 3 ) 2 = (pi + p 4 ) 2 = At 2 - 



(B.10) 
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This implies 

5^ = 18AV(/i 2 ) = - ln/i 2 + finite ) , (B.ll) 

where in the last equality we used dimensional regularization. To identify the beta function, 
consider the Callan-Symanzik equation 



G<"\{x,};ix,\)=0, (B.12) 



where fj, is the renormalization scale and G (jl \x\^ . . . , x n ) is the connected n-point function 
computed in renormalized perturbation theory. At one-loop order, 7(A) = 0. For the 
4-point function, 



4 

— 1 



G (i) {pi,P2,P3,P4:) = iM{pip 2 -> pm) IJ ^J, (B.13) 

we find 



9 „ m 27iA 2 A -i 



H^-G® = -^rW^- (B.14) 



Therefore the Callan-Symanzik equation is satisfied to order A 2 if 



q\2 

/5(A) = -^ + 0(A 3 ), (B.15) 

which is the familiar expression for the one-loop beta function. 

In the analogous computation for the renormalization of the double trace deformation, 
only the "s-channel" part of the above computation survives in the large N limit, as we 
have argued before. We write 

$ = $\*)+$\t) + $\u) (B.16) 
with 2 

5^\s) = 6A 2 VV) = -77^ (- - ln/x 2 + finite ) . (B.17) 

lbiT z \e ) 

Also at two-loop order, only "s-channel" diagrams will be important to leading order in 
1/N, as we have seen before; the relevant diagrams are in Figure [221 where the insertion of a 
cross denotes the one-loop counterterm 8^\s). We find (suppressing external propagators) 

= QtX + (6zA) 2 \iV(p 2 ) - iV(^ 2 )] - (6zA) 3 [V(p 2 ) - V(/2 2 )] 2 + . . . , (B.18) 

where p = p± + p%. The various V(fi 2 ) come from one-loop counterterms; a two-loop coun- 
terterm has merely cancelled a momentum-independent divergence (6zA) 3 [V(/i 2 )] 2 . Com- 
puting 

4 G ' 4)= ^-^ [l/(pVlV)1+ --- (ai9) 
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Figure 22: Two-loop diagrams in the large N limit. 



and 

= 6i - m\ [V(p 2 ) - IV)] + • • • , (B.20) 
it is easy to see that the Callan-Symanzik equation is satisfied for 

m = (B.21) 

It is easy to check that this remains the case at higher order in A. Thus we have verified 
that the beta function is one- loop exact. 

We will be interested in an approximation of the quantum effective action that is valid for 
a large range of field values, in particular for large field values. An appropriate framework 
is that of [50] , where the standard Feynman diagram expansion is resummed and the theory 
is organized in a derivative expansion, 

r = J d A x[-V{4>) - l -d^d^Z{$) + ...], (B.22) 

and an expansion in the number of loops. First, we focus on the effective potential, i.e., 
(minus) the term without derivatives, in the one-loop approximation. 
Starting from the classical potential 

Klass = ~\<P\ (B.23) 

one computes [50] 

^^d^vH 1 -^)' (a24) 

where the integral is over Euclidean four-momenta. This can be rewritten in terms of 
x = k 2 : 
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Since the integral diverges, we introduce a UV cutoff x < A 2 . For A > 0, we also need an IR 
cutoff x > n 2 IR if we want a real effective potential (see [64]). Thus we compute, dropping 
terms that vanish as A — > oo: 

A 2 / Q\^2\ Q\A2j,2 n\2j,4 ,,4 



x / 2 4 2 V ^jfl 

.^_, ?Ji) _^ ln f^!^).(B.26) 



tf R - 3A0 2 

For A < 0, i.e., a positive quartic potential, we can take fim = and recover the standard 
result [50] 

xr 1 / qXA2 , 2 9AV 9AV ^ -3A0 2 M 
^i-ioop = 7^ -3AA — + — — In — — . (B.27) 



32vr 2 [ ^4 2 V A 2 

For A > 0, the case of interest in the present paper, choosing fim = would give a complex 
effective potential, signaling the fact that localized wavefunctions are unstable to spreading 
while keeping the expectation value of <p fixed [64]. Therefore we choose instead to integrate 
out only the "unfrozen" (or "non-tachyonic" ) modes; thus we choose 

ti R = 3A</> 2 + e 2 (B.28) 

and let e — > at the end of the computation. This results in 

v 1 / <UA^2 9A 2 4 9A 2 4 /3A0 2 \1 

^i-ioop = |-3AA <P + — + — In ^— j j . (B.29) 

Now we add counterterms, imposing that the renormalized mass should vanish and defining 
the renormalized coupling A M by the renormalization condition 

V(Ai) = -^V- (B.30) 

This corresponds to choosing for the sliding scale \x the value of a spacetime-independent 
external scalar field <fi rather than, say, the momentum of an external line. This has the 
virtue that no large logarithms will interfere with perturbation theory, which will thus be 
reliable as long as the renormalized coupling A M is small. The result is 

A„ , 9A 2 4 ln(0 2 /u 2 ) 
V(<f>) = + " \ 4 l 2 ■ (B.31) 

The renormalization group equation can be obtained by demanding that V(<p) be indepen- 
dent of /j,: 

^ - "8^" (B ' 32) 

Here, we have ignored a contribution from d/dfi hitting the A 2 in the second term of (IB. 311) . 
which is justified as long as |A At ln(0///)| <C 1. Equation flB.32|) is solved by 



16tt 2 
9 liUjF/M 2 



^ = n,./..,/,^ . (B.33) 
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with M an arbitrary scale (this implements dimensional transmutation). Choosing fi = <fi, 
i.e., the renormalization scale is set by the value of the field 0, the Coleman- Weinberg 
potential can then be written as 



V(<f>) 



4tT 2 4 



(B.34) 



91n(0 2 /M 2 )' 



Now suppose that for some value O , the coupling is small, 



< < 1, 



(B.35) 



then M < \(j) \ and (1B.34I) is trustworthy (i.e., higher order corrections can be ignored) for 
any such that \<f>\ > |0o|. As a result, we can conclude that in this case 



(This also holds for the massive theory, as long asm< |0o|-) This analysis was for a single 
scalar field; as can be seen from the above discussion, it extends to the large N adjoint 
theory with only very minor changes. 

As we have seen, the one-loop diagram Figure [T8l involving two double trace vertices gives 
an order f 2 contribution to the beta function for the coupling /, leading to the effective 
potential ( 13.6j) . One can ask whether one-loop diagrams with one double trace vertex 
and one commutator squared vertex (see ( IB. II) ) could give contributions of order fg 2 to 
the beta function for / and thus to the effective potential, which would be important at 
small /. The answer is that the fg 2 contributions to the effective potential cancel. To 
see this, let us compute the one- loop effective potential for 0, where $ x (x) = <f>(x)U with 
TrU 2 = 1, in our deformed M = 4 SYM theory, see Q with (jElj) and (TO) . We choose 
a constant background $ x = cj)U and compute the masses M(0) of the various modes in 
this background; this involves fixing a gauge and introducing ghosts, see for instance [58]. 
Up to terms that can be absorbed in counterterms, every bosonic mode contributes 



to the effective potential, where A is a UV cutoff, while every fermionic mode contributes 
with a minus sign. In the undeformed M = 4 SYM theory, these contributions cancel 
exactly between bosons and fermions. Our double trace deformation changes the scalar 
masses but not the masses of gauge bosons and fermions, leading to a non-trivial effective 
potential. It is easy to see that the scalars § l ab have 



V(<P) 



oo for |0| 



oo. 



(B.36) 




(B.37) 




(B.38) 



{Ml h ) 2 = g 2 4> 2 {U aa -U bb ) 2 



2/oV 
N 2 



(l + 2S ab U 2 aa ). 



(B.39) 
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The contributions of order g A are those of the undeformed SYM theory and cancel by 
super symmetry against those of gauge fields and fermions. Contributions of order fg 2 
would have to come from off-diagonal scalars (a ^ b), but it is easy to see that they 
cancel between Q l and (i — 2, . . . 6), at least up to terms that can be absorbed in the 
counterterms. So we have verified that the one-loop contribution to the effective potential 
is at least of order f 2 . 

Since we will be interested in non-static solutions, knowledge of the potential alone is 
insufficient for our purposes: we also need to know the form of the derivative terms in 
the quantum effective action. Terms arising from non-divergent Feynman diagrams are 
suppressed by powers of the effective coupling A^, which is small in the regime where 
we trust our one-loop effective potential. It will turn out that such terms will be small 
for the configurations of interest. Logarithmically divergent Feynman diagrams induce, 
upon renormalization, terms involving ln(0//z), with \x the renormalization scale. Such 
terms disappear when we choose \i — <fi, as we have done for the effective potential. From 
power counting, the only possible divergences involving momenta arise in graphs with two 
external legs, which lead to wave function renormalization (as it turns out, such divergences 
do not occur at one- loop order in this theory, but do occur at higher loop order). We 
choose the renormalization condition Z(fx) = 1 in (IB. 221) . which for our choice p = 
implies that the field has a canonical kinetic term. Strictly speaking, the last conclusion 
holds for the theory of a single scalar field; in the adjoint model, the standard kinetic 
term — Tr9^$ (9^$ may be accompanied by terms with a different index structure, e.g. 
— (Tr$<9^$)(Tr$<9^$)/Tr$ 2 , which will, however, be suppressed by powers of the effective 
coupling (these terms do not arise form ultraviolet divergent Feynman diagrams, but from 
resumming diagrams with more than two external legs). 

The running coupling A^ given by (IB. 331) has a branch point at <fi = oo, the place where 
we want to impose boundary conditions that make the Hamiltonian self-adjoint. To avoid 
this complication, it is convenient to first study the system with a finite UV cutoff p on the 
momenta, in which case the potential will be meromorphic near = oo. The idea is that 
after imposing appropriate boundary conditions and computing particle creation, the UV 
cutoff can be removed and it can be checked that the results have a well-defined limit. So 
let us compute the one-loop effective potential, with all non-tachyonic ("unfrozen") modes 
integrated out, with a (smooth) UV cutoff p: 



l— loop 




(B.40) 



3A0V 




) 



i28vr 2 (9A 2 4 + p 4 ) 



Adding the tree-level potential and counterterms, we find 



V = --0 4 + Vx-ioop + A + B<j> 2 + C0 4 , 



(B.41) 
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where the last three terms are the counterterms. The coefficients A, B,C are fixed by the 
renormalization conditions 



V{0) = 0, V"(0) = 0, V{p) = -^-p\ (B.42) 



4 



where the dependence of the renormalized coupling A M on the renormalization scale p is 
indicated explicitly. Solving (IB. 421) . we find 

A = 0, (B.43) 
3A M p 2 



B 



1287T 



C = 128, W + «A^) h X >' - WW) - i P Hn (l - ) 

+ip 4 In (l + + 3?%^ ln(p 4 + 9A> 4 ) | . 

Substituting f)B.40|) and flB.431) in (1B.41|) . we obtain for the renormalized one-loop effective 
potential with renormalization scale p and UV cutoff p 

v = V | 3X ^ 2 p 2 

A V 128vr 2 p 2 (p 4 + 9A204) 
x {vrp 2 (p 4 + 9A> 4 ) - 9ttA> 4 2 - 6p 2 A,p 2 2 ln(3A,p 2 ) - zp 4 2 In (l - 

+ip A ( t> 2 In (l + ^-^j + 3p 2 2 A M p 2 ln(p 4 + 9Ajp 4 ) - vrp 2 p 4 

+6p%p 2 2 ln(3A M 2 ) + *p 4 p 2 In (l - - *p 4 p 2 In (l + 

-3pVA^ 2 ln(p 4 + 9Aj0 4 )}. (B.44) 

To complete the computation of the one-loop effective potential, we need A^ as a function 
of p, and we will eventually set p = <fi. For this reason, we impose that V should be 
p-independent at = p: 

= ^V\^. (B.45) 
This leads to the following equation for the running coupling at the scale p = 0: 

,dX(p 3p 2 A^ 



d<j> "" 16tt 2 2 (p 4 + 9A 2 4 ) 2 



x (B.46) 



x j-12p%0 2 + vrp 8 + 18vrp 4 Aj0 4 + 8l7rAj,0 8 - 54p 2 Aj0 6 ln (l + 
-*(p 8 + 27p 4 A 2 4 ) In (l - Jj£^ + *(p 8 + 27p 4 A 2 #) In (l + , 



P 4 



9XU 4 
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Figure 23: The running coupling A^ at finite UV cutoff p, versus z — 1/0, showing how it 
tends to the finite limit A^. 

C The Renormalized Effective Potential V{(p) in the Complex 0- 
plane 

The calculations performed in this paper rely on the use of complex solutions of the classical 
field equations for the effective theory. In the previous appendix, we computed the effective 
potential for real <fi by integrating out the non-tachyonic modes and renormalizing the 
theory at a scale defined by the value of the scalar field <fi. The effective potential is simply 
given by — ^A^ 4 , where A^, solves (1B.46I) with some given initial condition. One can also 
use (IB.46P to analytically continue V(<p) to complex <ft. If we take the p = oo limit before 
continuing, then V(<ft) has a branch point at = oo, which interferes with our construction 
of the self-adjoint extension via the method of images. Therefore we should rather continue 
at finite UV cutoff p, find the relevant classical solution and only then take the limit p — > oo. 

In order to study the form of the effective potential at large 0, it is convenient to set 
= 1/z. Equation flB.46h becomes 



dXff, 
dz 



3p 2 A z 2 



16tt 2 2 (A 4 + 9A2) 2 

-12p%;z 6 + irp 8 z 8 + 18vrp 4 A 2 ^ 4 + 81ttA 4 - 54p 2 A^ 2 In I 1 + 



(C.l) 



4 4 

p Z 



9Xt 



^(pV + 27p 4 A 2 z 4 )ln 1- 



ip 2 z 2 
3A^, 



+ *(pV + 27p 4 AV)ln( 1 + 



ip 2 z 2 ' 
3A,a 
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The right hand side is regular at z = 0, so given the coupling there, A^, it is straightforward 
to numerically integrate the equation to find A^ at real (f) (Figure [231) . Likewise, one can 
integrate in any particular direction in z to find the potential at complex <p. Since (1B.46[) 
is a first order differential equation, the solution A^ can only fail to be analytic at points 
where the right hand side is singular. Ignoring the singularity at small 0, where A^ becomes 
large and perturbation theory fails, the running coupling is singular only at solutions of 
p 4 + 9A^,0 4 = 0. Approximating A0 by its value at = oo, which is a good approximation 
when Aqo is small, one sees that these points are roughly located at (f> = p/y3X^, times a 
fourth root of unity. Clearly, as p tends to infinity, these four branch points converge to 
the origin, forming a single branch point at <f) = oo. 

The four singularities are all related, by two symmetries. First, the effective potential 
is real on the real z-axis, so by the Schwarz reflection principle its singularities occur at 
complex conjugate pairs. Second, as a result of the symmetry under cf) — ► —<p in the 
original (bare) theory, fIC.ll) is invariant under z —>■ — z. Thus we need only study a single 
branch point, for example the one at z — e ln ^ 4 '^3X^1 ' p = e l7r / 4 3s^/p, where we have defined 
3s| = A<£. 

Expanding the right hand side of ( 1C.1I) around z = e m ' 4 3s^/p, for example, we find 

- ^ \, ■ ,a — r4-^ + O Us* - e-^V/3)- 1 . (C.2) 

dz 128tt 2 (s - e-^pz/3) 2 vv H 1 ' > K ' 

At first sight, there appears to be a second order pole on the right hand side. But the 
situation is subtle because is itself a function of z. In the vicinity of the branch point, 
the solution is ~ s^^ + 5s(z), where z c = 3s ( f, )C e 1 '*/ 4 / p is the location of the branch point 
and s^c is the value of there. Setting \z\ = e x , Equation ( 1C.2|) becomes 



dSs _ 27s 5 c ln(2fa/e) 
~dx~ ~ 128n 2 5s 2 

which is solved for x approaching x c from below by 



(C.3) 



5s/ (He/(25sW « e-^(x c - xfs (j^j 3 . (C.4) 

Figure [23] shows the full numerical solution for s^, with the branch point being found by 
shooting, compared with this formula. We conclude that is finite at the branch points 
and varies in their vicinity roughly as Ss ~ (z c — z)s , with logarithmic corrections. The 
corresponding deviation of A^ from its value Aqo at infinite implied by (1Q.4[) is of order 

4 

A^ . Thus the fractional deviation in the coupling decreases as the running coupling to the 
one third power as p is raised. Therefore in the large p limit, the effective potential in the 
neighbourhood of each branch point becomes more and more accurately approximated by 
the effective coupling calculated at infinite p. 
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Figure 24: Behavior of the running coupling near a branch point in the complex z-plane. 
The third power of the left hand side of (1C.4I) is plotted against x c — x. The upper curve 
shows a linear fit. 

D Complex Classical Solutions in the Infinite Cutoff Limit 

In the previous appendix, we studied the behavior of the effective potential in the complex 

</>-plane for finite cutoff p. Besides the fourth order pole at <p = oo, the only singularities 

i 

are four mild, perturbatively suppressed branch points located at \z\b ~ A Jo /p. 

In Subsection 15.31 we explained how the coherent state parameter specifying the initial 
condition must be changed when a family of classical solutions crosses a branch point. 
We wish to follow the same procedure here, first at finite cutoff p and then taking the 
infinite p limit. As explained in Subsection 15.41 we start from classical solutions to the 
finite-p equations of motion which pass around the branch points without encountering 
any branch cuts so that the method of images may be used to construct the self-adjoint 
extension. These solutions have very large Xf, of order the cutoff in appropriate units, and 
times (tf — ti) which are roughly twice the classical rolling-down time. We now consider 
decreasing (tf — ti) so that the trajectory of the first classical solution, the upper solution 
in Figure is deformed towards the positive real x _ax is- As it does so, it encounters the 
two upper branch points and moves onto a new Riemann sheet. The associated shift in 
the coherent state parameter is calculated, as explained in Subsection 15.31 by comparing 
classical solutions passing to either side of the branch point. 

The question we want to address in this appendix is to what extent the shift in the co- 
herent state parameter can be calculated directly in the infinite-p limit. Because the branch 
points are very mild, it is plausible that the classical solutions of interest are insensitive 
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to the details of the finite-cutoff effective potential in their vicinity. Rather, the difference 
between neighbouring classical solutions which begin on different Riemann sheets accrues 
along the entire classical trajectory. But this latter effect can be computed by first taking 
p to infinity and then comparing the two classical solutions, one on the first Riemann sheet 
(on which the potential is real just above the positive real axis) and one on the second 
Riemann sheet (on which the potential is real just below the positive real axis). 

In this appendix we shall calculate the difference between two classical trajectories 
moving in slightly different potentials but with the same initial coherent state parameter 
(15. f p . the same final coordinate Xfi an d the same total time tf — tj. If we can show that 
the difference between these two trajectories is small, then we are justified in taking the 
infinite p limit before computing the shift in the coherent state parameter associated with 
moving onto the physical Riemann sheet. 

We shall work in the coordinate x used in Subsection 15.21 and the following subsections. 
We define the difference in potentials 5V(x) = V^x) — V p (x) where V p is the effective 
potential at UV cutoff p. We want to see whether the final classical solution is sensitive to 
the details of the potential at large values of x where the potential has significant dependence 
on p. In order to examine this sensitivity, we will solve the linearized equations for the 
perturbation to the classical solution induced by a change 5V(x) in the potential, for x 
near Xf. 

A small perturbation of the trajectory x(t) — > x(t) + 5x(t) which does not change the 
coherent state parameter, the final coordinate Xf and the total time tf — ti must obey the 
constraints 



_ 2 5e — V'(x)5xi 



5xi + 2iL 2 ^— i = 0, (D.5) 

hx; 



which relates 5xj to 5e, and 



St = 6 /" ^ = - f f d Je-V>(x)5x-5V(x) = 



where x = \/2{e — V{x)). 

The perturbation Sx(t) to the trajectory is found from the energy equation, 5e = x5x + 
V'(x)5x + SV(x) = 0. Imposing the final condition that 5xf = 0, this integrates to 

• f tf ~ 8V(x) . ft , 6e - SV(x) 
5x(t) = -x / dt v ; = -x / dx v ' . (D.7) 



Putting flD.51) and (ID. 71) together, we can calculate the perturbation in the initial position 
for the trajectory, 

2 /„, x ./-, x 



6x t = ^ I dx V -^ / dx - ." 3 v ; . (D.8) 



Exchanging orders in the integral, one finds 



<h , = it f Xf dx {Se - 5V(x))\n ((e - V{x))/{e - V{ Xi )))) ^ 
2 Li (2{e-V{x)))% 
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We are interested in the effect of a small change in the potential SV(x), in the region 
x ~ Xf affected by the cutoff, where the potential dominates over e. At large x the integral 
over 5V(x) looks like J dxx~ 2 lri2 (x)(SV/V) . But we have already argued, in the previous 
appendix, that for x in the vicinity of the branch points, SV/V is perturbatively small. 
It follows that the integral determining the shift in the initial condition 5xi is convergent 
at large x. In the limit as p goes to infinity, the fractional difference between the finite- 
and infinite-cutoff potentials, SV/V is localized at higher and higher values of x and its 
effect becomes negligible. Therefore the effect of moving to the next Riemann sheet may 
be computed using the infinite-cutoff effective potential. 
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